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PREFACE 


This  bulletin  contains  initial  results  and  analyses  derived  from  the  1983  MIZEX 
pilot  study  in  the  Greenland-Norwegian  Sea  in  June  and  July.  While  only  a  pilot 
study,  this  scientific  expedition,  carried  out  in  three  phases  with  two  personnel  ex¬ 
changes,  represented  one  of  the  most  complete  investigations  of  processes  at  a 
decaying  ice  edge  to  that  date.  Much  of  this  data  forms  an  independent  and  usable 
analysis  for  comparison  with  the  1984  main  summer  experiment.  This  report  lists  in¬ 
itial  results  from  several  of  the  experiments  with  special  emphasis  on  sea  ice 
dynamics  and  decay.  Other  analyses  will  appear  in  future  bulletins. 


4 


CONTENTS  ; 


{  Preface . 

Drift  Velocity  During  the  Drift-Station  Phase  of  MIZEX  83 

Miles  G.  McPhee. .  .... .... , , v, ,  _ _ _ 

^  MIZEX  83 — BIO  Buoy  Data  Summary  .  _ 

G.  Symonds and I.  Peterson.. . . 

|1  Analysis  ^ 


MIZEX  83  Mesoscale  Sea  Ice  Dynamics, 

W.D.  Hibler  III  and  M.  Lepparanta. . 

4  On  the  Rheology  of  a  Broken  Ice  Field  Due  to  Floe  Collision 
Hayley  Shen,  W.D.  Hibler  III  and  M.  Lepparanta _ _ 


<4Heat  and  Mass  Balance  Observations  During  the  MIZEX  83  Drift  Program 

Gary  A.  Maykut - _ - - - - - - - 

^  Ocean  Currents  and  Temperatures  in  the  Center  of  Fram  Strait  During 


MIZEX  83  j  —y 

Kenneth  Hunkins . .  . 

.Arctic  WhitecappingrPreliminary  Jffesults 
E.C.  Monahan  aril  P. A.  Bowyer . 


4— 


Page 

ii 

1 

13 

19 

29 


47 

53 


Drift  Velocity  During  the  Drift-Station  Phase  of  MIZEX  83 


MILES  G.  McPHEE 

McPhee  Research  Co.,  Rte  8,  Box  304H 
Yakima,  Washington  98908 


A  basic  initial  step  in  analyzing  measurements  of  velocity  taken  relative  to  a  drifting  ice 
station  is  to  determine  the  velocity  of  the  station  itself.  During  the  drift-station  phase  of 
MIZEX  83,  which  started  at  about  0900  on  day  179  and  ended  at  1200  on  day  189,  the  M/V 
PotarbjOm  was  moored  to  adjacent  ice  floes  and,  except  for  very  minor  adjustments,  its  tra¬ 
jectory  was  exactly  representative  of  those  floes.  The  geographic  location  of  the  ship  was 
monitored  by  satellite  positioning,  which  furnished  typically  about  40  fixes  per  day,  with  ac¬ 
curacy  somewhere  in  the  neighborhood  of  200  m. 

This  note  describes  a  fitting  technique  developed  for  on-site,  microprocessor  analysis  of 
station  velocity  to  be  implemented  for  the  1984  drift  experiment,  which  will  be  of  much 
longer  duration.  It  is  preliminary  in  the  sense  that  neither  a  rigorous  optimization  of  the  fit¬ 
ted  parameters  nor  a  velocity  error  analysis  has  been  carried  out.  Nevertheless,  it  presents  a 
potentially  useful  technique  that  has  not  appeared  elsewhere,  and  also  provides  a  velocity 
synthesis  that  may  be  helpful  for  investigators  who  had  fixed  and  profiling  current  measur¬ 
ing  systems  in  the  vicinity  of  the  ice  station.  The  appendix  contains  an  hourly  listing  of  the 
fitted  velocity  during  the  10  days  of  the  drift-station  phase. 

Suppose  that  over  some  period  of  time  we  can  approximate  an  ice  drift  trajectory  as  a  su¬ 
perposition  of  mean  motion  plus  oscillations  with  preferred  inertial  and  diurnal  frequencies. 
Differences  between  the  semidiurnal  and  inertial  frequencies  are  considered  unresolvable. 
The  oscillating  components  are  cyclical  (i.e.,  a  trajectory  due  to  the  oscillation  alone  returns 
to  its  initial  value  after  one  period),  but  not  necessarily  circular.  The  velocity  is  then; 

V(t)  -  Vm+S„t +D„e  ‘-’+Dmc‘~’  (1) 

where/is  the  angular  frequency  of  the  inertial  (semidiurnal  at  high  latitudes)  oscillation  with 
coefficients  S„  and  S^.,  the  subscripts  referring  respectively  to  clockwise  and  counterclock¬ 
wise  rotation.  fl„,  and  &>  similarly  describe  the  diurnal  oscillation. 

The  trajectory  for  an  arbitrary  time  reference,  t,  =  0,  with  initial  position  AT,,  is  given, 
from  integration  of  (1),  by 

*(/)  -  X,+  Vmt  +  (i/f)  15^e-«*-l)+5^1-e'^))+(//«)|f>^e  '--l)+De^l -e‘-*)l-  (2) 

For  illustration,  consider  a  trajectory  consisting  of  (a)  mean  motion  to  the  ENE  (60 0  bear¬ 
ing):  Ym  •  (8.66, 5)  cm  s'1;  (b)  a  pure  inertial  oscillation  with  amplitude  10  cm  s'1:  V,  » 
(-10,0)  cxp(-(ff);  and  (c)  a  diurnal  tidal  ellipse:  V,  -  (13,0)  exp(-/wf)  +  (0,5)  exp(iwr).  Note 
that  the  rotary  spectrum  of  such  motion  would  show  strong  peaks  in  the  clockwise  spectrum 
at /and  w,  with  a  weaker  peak  at  w  only  in  the  CCW  spectrum.  Figure  1  illustrates  the  two  os¬ 
cillating  trajectories  by  themselves,  and  then  the  total  trajectory  over  48  hr,  when  they  are  su¬ 
perimposed  on  the  steady  mean  drift.  Note  the  scale  changes  in  the  sketch. 


1.58  KM 


Fig.  1.  Trajectories  associated  with  a)  a  pure  inertial  oscillation 
with  period  12  hr  and  velocity  amplitude  10  cm/s;  (b)  a  tidal  ellipse 
with  period  24  hr  and  both  clockwise  and  counterclockwise  compo¬ 
nents;  and  fc)  a  superposition  of  (a)  and  (b)  with  a  mean  current  of 
10  cm/s  to  the  ENE. 


Determining  velocity  from  position  is  the  inverse  of  the  illustration  above.  The  problem 
may  be  stated  as:  from  position  estimates  at  irregularly  spaced  times,  estimate  the  six  coeffi¬ 
cients  that  give  the  best  fit  of  (2)  to  the  observations,  in  a  least-squares  sense,  over  some  time 
interval  comparable  to  the  period  of  the  lowest  frequency  oscillation.  The  technique  is  a  gen¬ 
eralization  of  one  introduced  by  Perkins  (1970)  in  the  analysis  of  inertial  waves  in  current 
meter  records  and  used  by  McPhee  (1980)  in  analyzing  AIDJEX  ice  drift.  In  essence,  it  is 
closely  related  to  estimating  spectral  density  at  discrete  frequencies  over  a  limited  time  peri- 
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od,  and  was  called  “complex  demodulation”  by  Perkins.  Naturally,  the  technique  works 
best  if  there  are  strong  signals  at  the  chosen  frequencies. 

Using  standard  least-squares  techniques  applied  to  complex  numbers,  the  problem  reduces 
to  solving  a  system  of  linear  equations. 

With  the  following  notation: 

B( a)  m  ( i/aXc—-l ) 

N 

E=  £ 

1-1 

where  N  is  the  number  of  position  samples  in  the  fitting  period,  the  equation  in  matrix  form 
is 
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LB{-u>)B{-  w)  D,„  LB[-u)X 


and  is  solved  by  Gaussian  elimination. 

Figure  2  shows  the  trajectory  of  the  M/V  Polarbjtrn  during  the  drift-station  phase  as  it 
was  plotted  on  the  ship.  Only  fixes  that  were  obviously  erroneous  were  omitted.  The  striking 
diurnal  loops  were  unexpected,  and  suggest  a  considerable  tidal  forcing,  since  the  inertial  pe¬ 
riod  at  81  °N  is  about  12.2  hr. 

Position  and  velocity  are  represented  in  a  cartesian  co-ordinate  system  consisting  of  a 
plane  centered  at  the  pole  with  x  increasing  along  a  line  projected  from  the  90  °W  meridian 
and  y  increasing  along  a  projection  of  the  Greenwich  Meridian.  Since  the  station  drift  was 
centered  near  7  °E  longitude,  positive  velocity  components  u  and  v  represent  approximately 
eastward  and  northward  motion  respectively.  The  transformation  from  geographic  co¬ 
ordinates  is 


X  =  (x  +  iy)  =  /?*(  sink  +  icosX) 

where  R  =  (90-lat)*111.2  km  and  X  is  the  longitude,  east  positive. 

The  fitting  procedure  described  in  the  previous  section  was  applied  to  the  data  shown  in 
Figure  2  in  two  passes.  First,  a  fit  for  four  of  the  components  (X0,  Vm,  S„,  and  S,cw)  was  cal¬ 
culated  every  3  hr  from  day  179.S  to  189.S  using  available  fixes  from  the  preceding  and  suc¬ 
ceeding  6  hr  (12-hr  data  window).  This  fit  passed  nearly  all  of  the  24-hr  energy  and  showed  a 
smaller  total  least-squared  error  than  a  similar  four-component  fit  for  the  diurnal  coeffi¬ 
cients.  The  original  SATNAV  fixes  were  then  compared  with  the  positions  as  calculated  by  the 
fit,  and  removed  from  the  data  time  series  whenever  they  deviated  by  more  than  500  m  from 


Fig.  2.  Plot  ofSATNA  V positions  as  received  at  the  Polarbj0rn  during  the  drift-station  phase 
of  the  MIZEX  83  experiment. 
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Fig.  3.  Fitted  trajectory  of  the  drift  station.  Each  dot  represents  1  hr,  with  pluses  at  half-daily 
intervals. 


the  fit.  On  the  second  pass,  the  data  (minus  the  out-of-bounds  points)  were  fitted  to  the  full  6 
(complex)-component  expression  using  the  inertial  period  (12.2  hr)  for  the  semidiurnal  com¬ 
ponents  and  24  hr  for  the  diurnal.  After  some  experimentation,  the  data  window  for  each  fit 
calculation  was  also  set  at  24  hr.  Fit  coefficients  were  calculated  every  3  hr.  All  the  calcula¬ 
tions  were  performed  on  an  Apple  II  microcomputer. 

Figure  3  shows  the  fitted  trajectory  with  time  marks.  Figure  4  shows  a  direct  comparison 
of  the  fitted  curves  with  all  observations  used  in  the  fits.  Figure  5  shows  (a)  the  total  velocity 
component,  where  the  heavy  line  represents  the  real  part,  u,  and  the  light  line  the  imaginary, 
v;  (b)  the  "mean”  velocity;  (c)  the  semidiurnal  part,  consisting  of 
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Fig.  4.  Comparison  of  fitted  positions  with  observations  used  in  calculating  the  fixes. 


V,  =  Sc„  exp(  -  ift)  +  S«.  exp(<//); 
and  (d)  the  diurnal  part, 


V,  -  D,w  exp(  -  iw  0  +  D«w  exp(/w  t). 

While  the  emphasis  in  this  note  is  not  on  interpretation,  the  semidiurnal  component  has 
been  labeled  “inertial,”  while  the  diurnal  is  called  “tidal.”  One  is  hard  pressed  to  explain  the 
24-hr  cycle  in  terms  other  than  tides,  and  the  lunar  phase  (full  on  25  June,  day  176;  last  quar¬ 
ter  on  3  July,  day  184)  is  consistent  with  the  attenuation  of  the  diurnal  signal  seen  in  the  first 
half  of  the  record.  The  semidiurnal  part  could  also  be  tidal,  but  seems  to  coincide  with  what 
small  wind  forcing  there  was  (evidenced  also  in  the  mean  velocity),  prompting  the  label  “in¬ 
ertial.”  Clearly,  that  identification  remains  tentative. 


Fig.  5.  (a)  Total  velocity,  u  lreal)-dark;  v  (imaginary)-Hght;  u  and  v  are  approximately 
eastward  and  northward  respectively;  (b)  mean  velocity,  V„;  (c)  semidiurnal  velocity 
(12. 2-hr  period);  and  (d)  diurnal  velocity. 
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APPENDIX 

Following  is  a  listing  of  hourly  values  of  position  (X)  in  km  and  velocity  ( K)  in  cm/s  in 
the  coordinate  system  described  in  the  text.  The  first  component  of  velocity  is  approxi¬ 
mately  eastward  positive;  the  second  is  approximately  northward. 
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(  20.88, 

12.18) 

Vt 

7 

( 

119.81,-971.61) 

(  17.76, 

10.12) 

8 

( 

120.40,-971.27) 

(  14.34, 

5.54) 

9 

c 

120.84,-971.14) 

(  10.58, 

1.76) 

10 

( 

121.17,-971.15) 

(  7.74, 

-1.84) 

11 

( 

121.35,-971.31) 

(  6.80, 

-4.72) 

12 

( 

121.60,-971.52) 

(  6.69, 

-6.82) 

■■  . 

13 

( 

121.84,-971.80) 

(  6.65, 

-8.41) 

14 

( 

122.12,-972.09) 

(  4.91,- 

-10.03) 

IS 

( 

122.26,-972.47) 

(  2.92,-10.66) 

16 

( 

122.33,-972.85) 

(  0.63,- 

-10.15) 

• 

17 

( 

122.31,-973.18) 

(  -1.71, 

-8.26) 

18 

( 

122.20,-973.43) 

(  -3.87, 

-5.61) 

19 

( 

122.03,-973.57) 

(  -5.54, 

-2.27) 

20 

( 

121.81,-973.60) 

(  -6.61, 

1.55) 

*  * 

21 

( 

121.57,-973.48) 

(  -6.49, 

5.25) 

22 

( 

121.36,-973.23) 

(  -5.31, 

8.40) 

23 

c 

121.19,-972.89) 

(  -3.19, 

10.54) 

24 

( 

121.13,-972.48) 

(  -0.32, 

11.99) 

\\V 

\V 

9 


I 

( 

121.17, 

-972.04) 

(  2.79, 

12.58) 

2 

( 

121.32, 

-971.59) 

(  5.45, 

12.51) 

3 

( 

121.56, 

-971.15) 

(  7.95, 

11.79) 

4 

( 

121.89, 

-970.74) 

(  9.80, 

10.61) 

5 

( 

122.26, 

-970.39) 

(  11.03, 

9.28) 

6 

( 

122.66, 

-970.09) 

(  11.33, 

7.58) 

7 

( 

123.06, 

-969.85) 

(  10.93, 

5.64) 

8 

( 

123.46, 

-969.69) 

(  10.02, 

3.15) 

9 

( 

123.80, 

-969.62) 

(  8.74, 

0.94) 

10 

( 

124.09, 

-969.62) 

(  7.49, 

-1.34) 

11 

( 

124.34, 

-969.70) 

(  6.39, 

-3.69) 

12 

( 

124.55, 

-969.88) 

(  5.52, 

-6.01) 

V 

13 

( 

124.74, 

-970.13) 

(  4.80, 

-8.04) 

14 

( 

124.90, 

-970.46) 

(  4.26, 

-9.80) 

15 

( 

125.03, 

-970.83) 

(  3.39, 

-10.35) 

16 

( 

125.13, 

-971.19) 

(  2.14, 

-9.83) 

■*/*  V- 

17 

( 

125.19, 

-971.52) 

(  0.69, 

-8.13) 

18 

( 

125.18, 

-971.77) 

(  -1.29, 

-5.60) 

19 

( 

125.10, 

-971.92) 

(  -3.28, 

-2.45) 

20 

( 

124.96, 

-971.95) 

(  -4.56, 

0.85) 

21 

( 

124.78, 

-971.85) 

(  -5.58, 

4.24) 

22 

( 

124.57, 

-971.64) 

(  -5.57, 

7.27) 

i 

23 

( 

124.38, 

-971.34) 

(  -4.38, 

9.74) 

24 

( 

124.27, 

-970.96) 

(  -1.73, 

11.37) 

Day  1 

L84 

1 

( 

124.27, 

-970.53) 

(  1.88, 

12.13) 

2 

( 

124.41, 

-970.10) 

(  5.75, 

12.39) 

3 

( 

124.69, 

-969.67) 

(  9.53, 

11.63) 

4 

( 

125.09, 

-969.27) 

(  12.54, 

10.26) 

_  . 

5 

( 

125.58, 

-968.94) 

(  14.65, 

9.09) 

6 

( 

126.12, 

-968.65) 

(  15.01, 

7.06) 

i 

7 

( 

126.65, 

-968.43) 

(  14.00, 

4.77) 

8 

( 

127.09, 

-968.31) 

(  11.38, 

2.40) 

9 

( 

127.46, 

-968.26) 

(  9.34, 

0.28) 

10 

( 

127.77, 

-968.29) 

(  7.49, 

-1.60) 

11 

( 

127.99, 

-968.37) 

(  5.33, 

-2.76) 

12 

( 

128.17, 

-968.49) 

(  4.52, 

-4.04) 

13 

( 

128.33, 

-968.66) 

(  4.16, 

-5.08) 

14 

( 

128.47, 

-968.86) 

(  3.86, 

-5.87) 

15 

( 

128.60, 

-969.08) 

(  3.01, 

-6.07) 

16 

( 

128.68. 

-969.29) 

(  1.36, 

-5.65) 

17 

( 

128.67, 

-969.47) 

(  "2.46, 

-4.71) 

18 

( 

128.53, 

-969.60) 

(  -5.06, 

-2.91) 

19 

( 

128.31, 

-969.67) 

(  -7.43, 

-0.40) 

20 

( 

128.00, 

-969.62) 

(-10.45, 

2.71) 

21 

( 

127.61, 

-969.47) 

(-11.03, 

5.88) 

22 

( 

127.22, 

-969.20) 

(-10.10, 

8.84) 

1 

23 

( 

126.90, 

-968.84) 

(  -9.13, 

11.54) 

24 

( 

126.63, 

-968.40) 

(  -5.75, 

13.12) 

Day  185: 

1 

( 

126.49, 

-967.91)  i 

(  -1.72, 

13.65) 

2 

( 

126.48, 

-967.43) 

(  1.03, 

13.32) 

3 

( 

126.59, 

-966.97) 

(  4.26, 

12.04) 

4 

( 

126.78, 

-966.57) 

(  6.27, 

10.16) 

5 

( 

127.00, 

-966.25) 

(  5.17, 

8.27) 

6 

( 

127.18, 

-965.99) 

(  4.29, 

6.18) 

i 

7 

( 

127.30, 

-965.80) 

(  2.42, 

4.23) 

8 

( 

127.36, 

-965.68) 

(  -0.20, 

2.54) 

9 

( 

127.29, 

-965.61) 

(  -3.23, 

1.10) 

10 

( 

127.13, 

-965.59) 

(  -6.13, 

-0.07) 

11 

( 

126.82, 

-965.61) 

(  -9.44, 

-1.32) 

12 

( 

126.46, 

-965.67) 

(-10.57, 

-2.09) 

.  *  • 

13 

( 

126.07, 

-965.76) 

(-10.79, 

-2.62) 

14 

( 

125.68, 

-965.85) 

(-11.55, 

-3.34) 

•  \  *  ^ 

15 

( 

125.28, 

-965.97) 

(-10.90, 

-3.41) 

16 

( 

124.89, 

-966.09) 

(-10.37, 

-2.99) 

17 

( 

124.50, 

-966.18) 

(-11.88, 

-2.35) 

18 

( 

124.06, 

-966.24) 

(-12.25, 

-0.86) 

19 

( 

123.61, 

-966.24) 

(-12.90, 

1.14) 

20 

( 

123.14, 

-966.16) 

(-14.23, 

3.41) 

21 

( 

122.61, 

-965.99) 

(-15.06, 

5.76) 

22 

( 

122.06, 

-965.75) 

(-15.52, 

7.82) 

23 

( 

121.50, 

-965.42) 

(-15.41, 

10.09) 

24 

( 

120.95, 

-965.04) 

(-14.66, 

10.62) 

Day  186: 

i 

( 

120.45, 

-964.67) 

(-13.28, 

10.16) 

2 

( 

120.00, 

-964.31) 

(-10.06, 

9.54) 

3 

( 

119.67, 

-964.00) 

(  -8.43, 

7.67) 

4 

( 

119.39, 

-963.76)  i 

(  -7.36, 

5.67) 

■  ■  . 

5 

( 

119.17, 

-963.57) 

(  -5.77, 

4.09) 

6 

( 

118.96, 

-963.45) 

(  -6.07, 

3.10) 

7 

( 

118.73, 

-963.34) 

(  -6.68, 

3.01) 

8 

( 

118.49, 

-963.22)  i 

!  -5.94, 

3.45) 

9 

( 

118.27, 

-963.07) 

(  -6.11, 

4.68) 

10 

( 

118.05, 

-962.88) 

(  -5.75, 

6.00) 

11 

( 

117.89, 

-962.64) 

(  -3.22, 

7.75) 

12 

( 

117.80, 

-962.36)  i 

(  -1.93, 

7.62) 

13 

( 

117.75, 

-962.10) 

(  -0.61, 

6.20) 

14 

( 

117.78, 

-961.92) 

(  2.83, 

4.23) 

15 

( 

117.88, 

-961.82) 

(  2.70, 

1.19) 

16 

( 

117.96, 

-961.83) 

(  1.43, 

-1.67) 

17 

( 

118.02, 

-961.93) 

(  1.37, 

-3.35) 

18 

( 

118.02, 

-962.06) 

(  -1.24, 

-3.57) 

19 

( 

117.93, 

-962.16) 

(  "3.64, 

-1.89) 

20 

( 

117.80, 

-962.17)  i 

(  -3.9' , 

1.21) 

» 

21 

( 

117.65, 

-962.05) 

(  -4.08, 

5.88) 

22 

( 

117.52, 

-961.74) 

(  -2.53, 

10.91) 

23 

( 

117.49, 

-961.28) 

(  0.66, 

15.11) 

24 

( 

117.58, 

-960.68) 

(  4.91, 

18.03) 

MIZEX  83— BIO  Buoy  Data  Summary* 


G.  SYMONDS  AND  I.  PETERSON 

Bedford  Institute  of  Oceanography,  Atlantic  Oceanographic  Laboratory 
Dartmouth,  Nova  Scotia,  Canada 


This  preliminary  report  contains  a  summary  of 
the  trajectories  and  velocities  of  the  six  drifting 
buoys  deployed  by  the  Bedford  Institute  of  Ocean¬ 
ography  during  MIZEX  83.  The  buoys  were 
tracked  using  the  ARGOS  positioning  system. 
Dates  and  positions  of  deployment  and  recovery 
are  listed  in  Table  1.  Buoy  2404  ceased  operating 
shortly  after  deployment  and  is  not  included  in  the 
following  figures.  Buoys  2405  and  2406  were  not 
recovered.  Buoy  2406  ceased  operating  on  24  July 
1983,  while  buoy  2405  continued  until  19  August 
1983.  Buoys  2407  and  2409  were  recovered  at  the 
end  of  the  drift  phase  and  redeployed  in  the  East 
Greenland  Current  during  the  edge  phase  (14-23 
July  1983).  Buoy  2408  was  recovered  by  R/V 
Lance  on  21  July  1983. 

Also  listed  in  Table  1  are  the  standard  devia¬ 
tions  in  the  u  (eastward)  and  v  (northward)  com¬ 


ponents  of  velocity.  In  this  report,  velocity  is  com¬ 
puted  from  1  -hr  linearly  interpolated  positions, 
and  then  smoothed  with  a  2-hr  running  mean.  The 
relatively  high  standard  deviation  for  buoy  2409  is 
due  largely  to  instrument  error,  since,  at  a  fixed 
location  on  land  before  the  field  experiment,  the 
standard  deviation  of  position  for  buoy  2409  was 
about  twice  that  for  the  other  buoys. 

Mean  trajectories  were  computed  from  2-hr  lin¬ 
early  interpolated  positions  with  a  24-hr  running 
mean.  Trajectories  of  the  three  buoys  deployed 
prior  to  the  drift  phase  (2405,  2406,  and  2408)  are 
plotted  in  Figure  1 .  The  trajectories  of  buoys  2407 
and  2409  are  shown  both  during  the  drift  phase 
(Figure  2)  and  later  during  the  edge  phase  (Figure 
3). 

A  listing  of  24-hr-averaged  positions  for  each 
buoy  is  given  in  Table  2. 


Table  1.  Deployment  and  recovery  dates  and  positions  and  standard  deviations  of  velocity 
for  the  BIO  buoys. 


Buoy 

no. 

Deployment 

Recovery 

Length 

(days) 

S.D.  temts) 

Date 

Position 

Date 

Position 

u 

V 

2404 

18  Jun  83 

80.68  *N,  12.09  "E 

20  Jun  83 

40.5114,  12.69  °E 

3 

2405 

18  lun  83 

81.37  °N,  10.60  “E 

19  Aug  83 

80.6914,  5.92  °E 

63 

09.38 

7.78 

2406 

18  Jun  83 

81.0214, 09.08  ^ 

24  Jul  83 

80.5214,  11.03'E 

37 

9.87 

8.82 

2407 

18  Jun  83 

80.6014,  11.61  "E 

21  Jun  83 

80.33 °N,  ll.76°E 

4 

— 

— 

2407 

27  Jun  83 

80.86  14,  08.551: 

10  Jul  83 

81.0314,  V).i2"t 

14 

8.66 

8.30 

2407 

14  Jul  83 

79.3514,  0I.44TE 

22  Jul  83 

78.29 “N,  03.801; 

9 

— 

— 

2406 

18  Jun  83 

80.8314,  13.00  "E 

21  Jul  83 

80.5314,  21.11  °E 

34 

12.79 

10.76 

2409 

27  Jun  83 

81.0614,  08.62  "E 

8  Jul  83 

81.0914,  09.80°E 

12 

13.47 

9.89 

2409 

14  Jul  83 

79.36  14,  0I.82°E 

22  Jul  83 

78.5719,  03.50  °E 

9 

_ 

— 

•This  report  is  meant  only  as  a  very  brief  summary  for  inter¬ 
ested  MIZEX  (dentists  and  should  not  be  considered  a  final 
scientific  publication. 
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MIZEX  83  Mesoscale  Sea  Ice  Dynamics:  Initial  Analysis 


W.D.  HIBLER  III 

U.S.  Army  Cold  Regions  Research  and  Engineering  Laboratory 
Hanover,  New  Hampshire  03755-1290 

M.  LEPPARANTA 
Institute  of  Marine  Research 
Helsinki,  Finland 


This  report  describes  the  basic  characteristics  of 
the  mesoscale  sea  ice  dynamics  data  set  taken  dur¬ 
ing  the  drifting  phase  of  the  MIZEX  83  pilot  study. 
The  basic  data  set  consists  of  distance  measure¬ 
ments  accurate  to  1-2  m  between  Del  Norte  units 
at  five  sites  with  one  site  being  the  ship.  In  addi¬ 
tion  to  these  measurements,  the  overall  rotation  of 
the  array  was  determined  with  less  accuracy  by  us¬ 
ing  the  ship's  radar  to  measure  the  angles  of  either 
active  or  passive  radar  transponders  on  the  ice  at 
2-hr  intervals. 

The  basic  configuration  of  the  array  is  shown  in 
Figure  1,  which  is  drawn  to  scale  for  29  June  0800 
GMT.  Each  complete  Del  Norte  measurement  con¬ 
sists  of  measuring  directly  the  7  distances  d„,  d„ 
du  dit  b„  bi,  bi  as  noted  in  this  figure.  These  dis¬ 
tances  are  measured  to  an  acccuracy  of  1-2  m. 
Other  distances  are  triangulated  and  will  be  less 
accurate,  with  the  precise  measurement  error  de¬ 
pending  on  the  geometry  of  the  array.  It  was  ini¬ 
tially  planned  to  make  measurements  every  hour 
or  so,  but  once  data  began  to  be  acquired  it  be¬ 
came  obvious  that  shorter  time  intervals  would  be 
necessary  to  resolve  all  motion  occurring.  Accord¬ 
ingly,  the  acquisition  computer  was  repro¬ 
grammed  to  automatically  acquire  data  at  3-min 
intervals.  This  operational  mode  was  employed 
throughout  the  experiment  although  there  were 
occasional  data  gaps  due  to  battery  failure,  polar 
bears,  and  interference  from  less  accurate  Motor¬ 
ola  transponders  located  on  certain  floes.  Such  ac¬ 
tive  “ship’s  radar”  transponders  will  not  be 
necessary  in  the  future,  as  it  was  found  that  pas¬ 
sive  transponders  could  be  effectively  located  for 
angular  measurements  provided  their  distances 
from  the  Del  Norte  units  were  known.  It  should 
also  be  noted  that  the  rotation  of  the  array  is  criti¬ 
cal  only  when  one  wants  the  north  and  east  com- 


i _ i  i  i _ i  i 

0  2  4  6  8  (Ohm 


Fig.  1.  Configuration  of  strain  array 
and  nomenclature  for  distances  dis¬ 
cussed  later  in  the  paper.  The  triangle 
is  plotted  to  scale  for  29  June  0800 
GMT;  the  ship  (M/V  Polarbj  firn)  is 
located  at  the  center  point  P  moored  to 
an  ice  floe.  Del  Norte  units  are  located  at 
circled  points.  Directly  measured 
distances  (accurate  to  1-2  m)  are  denot¬ 
ed  by  subscripted  b's  and  d's  while  dis¬ 
tances  calculated  by  triangulation  (less 
accurate  depending  on  geometry)  are  de¬ 
noted  by  subscripted  s 's. 

ponents  of  relative  ice  drift.  The  Del  Norte  system 
itself  gives  the  strain  rate  invariants  very  accu¬ 
rately  without  the  need  for  any  other  angular 
measurements. 

While  the  main  measurement  program  occurred 
with  the  ship  moored  to  a  floe,  data  were  also 
taken  before  and  after  the  center  floe  deployment 
with  the  ship  moving.  This  was  made  possible  by 
the  automated  nature  of  the  system,  which  al¬ 
lowed  a  measurement  to  be  made  in  about  IS 
seconds,  during  which  time  the  ship’s  motion  was 
relatively  small.  T riangulation  could  then  be  made 
of  the  relative  positions  of  points  on  the  ice. 


Fig.  2.  Photograph  of  ice  floes  near  the  center  of  the  strain  array.  The  location  of  the  Polar- 
bj0rn  is  noted,  which,  for  scale,  is  50  m  long.  Photograph  taken  by  Dr.  Vernon  Squires  of 
the  Scott  Polar  Research  Institute  as  part  of  an  aerial  survey  carried  out  by  helicopter  flights. 


In  terms  of  the  overall  field  operation,  there 
were  no  major  problems,  and  generally  the  system 
operated  much  better  than  expected.  Battery  life¬ 
times  were  about  4  days,  but  varied  with  different 
types  of  transponders.  The  connection  system  for 
the  batteries  will  be  improved  for  MIZEX  84  to  fa¬ 
cilitate  battery  replacement.  One  minor  problem 
was  that  the  tripod  height  (about  4  m)  was  a  bit 
low  for  the  initial  baseline  unit,  which  was  initially 
at  site  C.  To  correct  this,  the  baseline  tripod  was 
raised  about  1.5  m  and  moved  to  site  D.  After 
about  5  days  one  master  unit  began  malfunction¬ 
ing  and  had  to  be  replaced,  but  overall  the  auto¬ 
mated  nature  of  the  system  together  with  its  high 
accuracy  proved  most  valuable.  These  features  al¬ 
lowed  the  acquisition  of  by  far  the  most  detailed 
and  temporally  dense  ice  deformation  measure¬ 
ments  in  existence  to  date. 


In  addition  to  transponders,  sites  O,  C  and  J 
were  equipped  with  current  meter  strings  moored 
to  toroid  buoys.  As  discussed  below,  the  large  de¬ 
formation  rates  observed  demonstrate  the  neces¬ 
sity  of  accurate  ice  deformation  measurements  in 
order  to  precisely  interpret  the  current  measure¬ 
ments.  Basically,  relative  ice  drift  rates  of  2  to  3 
cm/s,  and  sometimes  8  cm/s,  were  common. 
Often  .hese  events  would  last  for  only  a  few  hours, 
so  fine  temporal  resolution  becomes  critical.  Be¬ 
cause  of  these  large  ice  motion  fluctuations,  ob¬ 
served  variations  in  current  meter  readings  may  in 
many  cases  be  due  to  ice  motion. 

The  ice  pack  was  relatively  homogeneous  over 
the  region  covered  by  the  triangle.  Ice  floes  were 
typically  30-50  m  in  diameter,  and  the  compact¬ 
ness  of  the  ice  was  80  to  90% .  Figure  2  shows  a 
photograph  over  the  ice  taken  from  the  helicopter 


a 


180  182  184  186  188  190 

Julion  Doy 

Fig.  3.  The  area  of  the  triangle  OCJ  at  ‘A -hr  intervals 
filtered  with  a  tow-pass  filter  with  transition  band  from 
48-  to  60-minute  periods.  The  high-frequency  variations 
during  days  184.5  to  186.5  are  due  to  problems  with  the 
measurement  system  (a  poorly  functioning  master  unit). 


at  an  altitude  of  about  5000  ft.  To  a  first  approxi¬ 
mation  the  ice  pack  is  an  almost  ideal  two-dimen¬ 
sional  granular  medium,  and  our  measurement 
scale  is  good  for  studying  its  small-scale  deforma¬ 
tion.  The  pack  was  a  mixture  of  first-year  and 
multiyear  ice,  with  thickness  ranging  from  1  to  4 
m. 

The  overall  long-term  character  of  the  deforma¬ 
tion  is  shown  in  Figures  3  through  5.  Figure  3 
shows  the  area  of  the  triangle  at  '/i-hr  intervals 
over  the  drift  station  experiment.  The  curve  has 
been  smoothed  with  a  low-pass  symmetric  filter 
with  241  weights  and  a  transition  band  from  48-  to 
60-min  periods.  The  filter  design  is  described  in 
Hibler  (1972).  Figure  4  shows  the  triangle  config¬ 
uration  at  approximately  the  beginning,  middle, 
and  end.  Further,  Figure  4  shows  the  approximate 
relative  position  and  orientation  of  the  ice  edge. 
As  can  be  seen  from  Figure  3,  the  area  shows  a 
rapid  opening  and  closing  early  in  the  experiment 
followed  by  a  gradual  convergence  and  then  a  ma¬ 
jor  closing-opening-closing  cycle  during  the  last 
several  days  of  the  experiment.  Examination  of 
Figure  4  shows  that  this  event  over  the  last  several 
days  is  accompanied  by  a  major  shearing  of  the 
triangle,  with  the  southern  points  moving  north¬ 
eastward  relative  to  the  northern  point.  This  over¬ 
all  convergence  is  consistent  with  the  light  south¬ 
erly  winds  during  the  experiment  and  the  generally 
northeastward  drift,  which  would  tend  to  push  the 
ice  edge  back  into  the  main  pack.  Note,  however, 
that  there  are  fluctuations  on  top  of  this  trend. 


0  2  4  6  8  10km 
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- 8  Jul  (0800  GMT) 


Open  Woter 

Fig.  4.  Configuration  change  of  the  triangle 
over  the  period  of  the  experiment.  Also 
shown  is  the  approximate  location  and  ori¬ 
entation  of  the  ice  edge  drawn  to  scale. 


Moreover,  although  not  shown  there  is  some  indi¬ 
cation  that  the  shearing  is  greater  than  the  con¬ 
vergence,  as  one  would  expect  for  a  relatively  in¬ 
compressible  two-dimensional  interacting  ice 
field. 
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Fig.  5.  Orientation  of  the  baseline  (line  from  the  ship  to 
site  D)  and  the  heading  of  M/  V  Polarbj  0rn  at  I -hr  inter¬ 
vals  over  the  period  of  the  experiment.  The  angle  in¬ 
creases  counterclockwise,  and  the  zero  points  east. 


With  regard  to  orientation,  the  baseline  (M/V 
Polarbj  dm — site  D)  rotated  about  20 0  during  the 
drifting  phase  (Fig.  5).  This  caused  a  relative  mo¬ 
tion  in  the  array  comparable  to  the  strain  motion 
measured  with  the  Del  Norte  system.  Because  the 
accuracy  of  the  orientation  angle  is  only  0.5 
large  uncertainties  arise  when  transforming  the 
Del  Norte  data  to  a  north-east  coordinate  system. 
These  uncertainties  do  not,  however,  affect  the 
strain  rate  invariants.  At  the  distance  of  5  km,  a 
0.5  °  sector  is  40  m  wide.  Also  shown  in  Figure  5  is 
the  rotation  of  the  Polar bj dr n  measured  with  a 
gyrocompass.  The  large  bump  in  day  184  is  man- 
caused  because  of  the  necessity  to  change  the 
ship’s  mooring.  The  two  rotation  curves  show 
similar  behavior,  but  the  high-frequency  fluctua¬ 
tions  are  much  more  pronounced  for  Polarbj drn. 

Table  1  gives  the  position  of  the  Del  Norte  sites 
in  north-east  coordinates  at  6-hr  intervals.  Their 
accuracies  are  6x  s  0.01  •  \y\,  6y  =0.01  •  |x|. 
These  large  errors  are  due  to  the  orientation  uncer¬ 
tainty  discussed  above.  The  baseline  orientation  is 
given  in  Table  2.  Note  that  the  direction  is  counted 
counterclockwise  from  the  east;  the  numbers  thus 
tell  that  the  direction  is  slightly  southward  from 
west  (see  also  Fig.  1).  The  data  in  Table  1  were  used 
to  estimate  the  zeroth  and  first-order  terms  of 
6-hour  mean  ice  velocities  using  the  linear  model 
(see  e.g.  Hibler  et  al.,  1974) 

u  -  u0+x*d,u|o  +  ydu\*  +  e., 

v  =  v,  +  x*d,v|0  +  yd  v j0  +  e„ 


where  the  subscript  0  refers  to  Polarbjdrn  and  (e„, 
ev)  is  the  error  term.  The  results  are  presented  in 
Table  3.  The  inaccuracy  due  to  measurement  error 
is  0. 1  cm/s  for  the  constant  terms  u0  and  v0  and 
0.1%  h'1  for  the  first-order  derivatives.  Also 
shown  are  the  principal  strain  rates  (El  and  E2), 
the  direction  of  the  first  principal  axis  (AE1),  di¬ 
vergence  (El),  twice  the  maximum  shear  rate 
(Eli),  and  vorticity  (VOR).  The  goodness  of  fit, 
defined  as  the  ratio  of  the  explained  variance  to 
the  total  variance  (last  column),  is  typically  85- 
95%. 

Analysis  of  finer  time-scale  data  shows  that  a 
large  variety  of  significant  fluctuations  is  superim¬ 
posed  on  the  overall  trends.  This  is  illustrated  in 
Figures  6  and  7,  which  show  raw  and  filtered  dis¬ 
tances  at  3-min  intervals.  As  can  be  seen,  the  fil¬ 
tered  record  shows  large  variations  on  time  scales 
of  a  few  hours.  In  the  raw  data,  there  are  signifi¬ 
cant  fluctuations  over  ‘A -hr  intervals.  These  fluc¬ 
tuations  are  undoubtedly  related  to  floe  bumping 
as  well  as  randomness  in  the  forcing  fields.  How¬ 
ever,  it  is  also  possible  that  they  are  due  to  more 
organized  floe  interactions,  which  can  manifest 
themselves  as  kinematic  waves  propagating  across 
the  array.  Numerical  studies  by  Hibler  et  al.  (1983) 
have  shown  that  such  waves  and  fluctuating  ef¬ 
fects  in  time  can  arise  in  a  nonlinear  plastic  sea  ice 
model  even  though  the  external  forcing  is  tempor¬ 
ally  constant.  The  simplest  example  of  such  a 
wave  is  a  "ridging  front”  that  can  occur,  for  ex¬ 
ample,  if  one  pushes  in  one  end  of  a  series  of 
equally  separated  floes  (or  billiard  balls);  the  colli- 
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Fig.  6.  Raw  measured  distances  during  a  12-hr  period  in  the  early  part  of 
the  experiment.  The  vertical  scale  shows  the  distance  scale,  and  absolute 
initial  values  are  listed  for  each  curve. 


i- 


.  f. 


Fig.  7.  Distances  filtered  with  a  low-pass  filter  (transition  band  48-  to 
60-minute  periods)  during  a  30-hr  period  in  the  late  part  of  the  experiment. 
The  vertical  scale  shows  the  distance  scale,  and  absolute  initial  values  are 
listed  for  each  curve. 
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Tabic  1.  Positions  of  the  Del  Norte  stations  with  respect  to  Polarbjfrn  in  north-east  coordinate  sys¬ 
tem.  x  and  y  are  kilometers. 


Station 


Cyclesonde 


1  Mo 

Day 

Hr 

(GMT) 

X 

y 

X 

y 

X 

y 

X 

y 

:  6 

28 

1800 

-1.704 

-0.153 

-2.701 

-4.586 

-2.164 

5.237 

4.790 

-0.584 

6 

29 

000 

-1.616 

-0.213 

-3.124 

-3.957 

-2.217 

4.645 

5.417 

0.661 

6 

29 

600 

-1.702 

-0.439 

-3.206 

-4.110 

-2.894 

3.982 

4.816 

-0.511 

—  - - 

1  6 

29 

1200 

-1.655 

-0.473 

-3.105 

-4.212 

-2.835 

3.953 

4.682 

-0.405 

•  , 

6 

29 

1800 

-1.676 

-0.508 

-3.046 

-4.216 

-3.139 

3.766 

4.820 

-0.337 

6 

30 

000 

-1.717 

-0.408 

-3.126 

-3.956 

-3.364 

3.735 

4.884 

-0.883 

6 

30 

600 

-1.714 

-0.389 

-3.121 

-3.932 

-3.329 

3.777 

4.848 

-0.968 

6 

30 

1200 

-1.674 

-0.386 

-3.100 

-3.962 

-2.997 

3.857 

4.727 

-0.834 

6 

30 

1800 

-1.683 

-0.338 

-3.051 

-3.866 

-3.229 

3.936 

4.783 

-1.054 

7 

1 

000 

-1.673 

-0.350 

-3.091 

-3.874 

-3.096 

3.907 

4.743 

-0.955 

3  7 

1 

600 

-1.653 

-0.317 

-3.148 

-3.846 

-3.006 

3.992 

4.607 

-1.101 

• 

1  7 

1 

1200 

-1.626 

-0.306 

-3.145 

-3.840 

-2.882 

4.059 

4.590 

-1.098 

7 

1 

1800 

-1.631 

-0.294 

-3.204 

-3.803 

-2.876 

4.106 

4.399 

-1.490 

:•  7 

2 

000 

-1.631 

-0.322 

-3.047 

-3.887 

-3.016 

4.123 

4.483 

-1.232 

-  7 

2 

600 

-1.638 

-0.304 

-3.108 

-3.820 

-2.932 

4.097 

4.527 

-1.172 

7 

2 

1200 

-1.628 

-0.305 

-3.104 

-3.825 

-2.907 

4.105 

4.492 

-1.163 

7 

2 

1800 

-1.637 

-0.281 

-2.985 

-3.823 

-3.076 

4.087 

4.419 

-1.747 

" 

i  7 

3 

000 

-1.628 

-0.336 

-2.887 

-3.903 

-3.128 

3.999 

4.484 

-1.455 

s _ 

1  7 

3 

600 

-1.655 

-0.292 

-2.986 

-3.805 

-3.071 

4.022 

4.554 

-1.483 

7 

3 

1200 

-1.649 

-0.287 

-2.974 

-3.804 

-3.095 

4.018 

4.473 

-1.638 

v'v 

7 

3 

1800 

-1.679 

-0.279 

-2.929 

-3.791 

-3.218 

4.048 

4.555 

-1.746 

.•.V.*  \ 

7 

4 

000 

-1.663 

-0.315 

-2.911 

-3.831 

-3.198 

3.965 

4.418 

-1.967 

7 

4 

600 

-1.688 

-0.314 

-2.984 

-3.814 

-3.153 

3.904 

4.400 

-1.840 

7 

4 

1200 

-1.664 

-0.396 

-3.057 

-3.914 

-3.175 

3.788 

4.240 

-1.632 

*  • 

■  7 

4 

1800 

-1.654 

-0.335 

-3.211 

-3.837 

-3.016 

3.919 

4.199 

-1.586 

1  7 

5 

000 

-1.638 

-0.463 

-2.975 

-4.085 

-3.178 

3.685 

4.184 

-1.541 

7 

5 

600 

-1.650 

-0.432 

-3.239 

-3.931 

-3.091 

3.696 

4.048 

-1.891 

7 

5 

1200 

-1.655 

-0.386 

-2.936 

-4.041 

-3.129 

3.755 

4.115 

-1.846 

7 

5 

1800 

-1.595 

-0.475 

-2.675 

-4.163 

-3.459 

3.632 

4.318 

-1.186 

7 

6 

000 

-1.572 

-0.521 

-2.431 

-4.224 

-3.725 

3.474 

4.448 

-0.993 

7 

6 

600 

-1.576 

-0.545 

-2.367 

-4.256 

-3.800 

3.375 

5.283 

-1.005 

,  *  ;  .  ' 

|  7 

6 

1200 

-1.671 

-0.476 

-2.325 

-3.879 

-4.152 

3.361 

5.276 

-1.407 

• 

1  7 

6 

1800 

-1.715 

-0.626 

-2.208 

-3.947 

-4.306 

2.966 

5.507 

-1.683 

7 

7 

000 

-1.683 

-0.641 

-2.149 

-3.959 

-4.300 

2.900 

5.393 

-1.508 

7 

7 

600 

-1.713 

-0.742 

-2.075 

-3.975 

-4.458 

2.666 

5.817 

-1.484 

•’v  V  ' 

7 

7 

1200 

-1.694 

-0.808 

-2.054 

-4.034 

-4.442 

2.529 

5.925 

-1.097 

7 

7 

1800 

-1.686 

-0.647 

-2.084 

-3.773 

-4.653 

2.787 

5.872 

-1.738 

7 

8 

000 

-1.598 

-0.548 

-2.267 

-3.798 

-5.094 

3.142 

5.329 

-2.074 

*.‘v  ’ 

7 

8 

600 

-1.443 

-0.389 

-1.996 

-3.420 

-5.210 

3.180 

5.365 

-1.907 

• 
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Table  2.  Orientation  of  the  baseline  in 
north-east  coordinate  system.  The  angle  is 
counted  counterclockwise  from  the  east. 


Dace 

Angle 

Mo 

Day 

Hr 

(GMT  ) 

6 

28 

1200 

189,9 

6 

28 

1800 

185.1 

6 

29 

000 

187.4 

6 

29 

600 

194.5 

6 

29 

1200 

195.9 

6 

29 

1800 

196.9 

6 

30 

000 

193.4 

6 

30 

600 

192.8 

6 

30 

1200 

193.0 

6 

30 

1800 

191.4 

7 

1 

000 

191.8 

7 

1 

600 

190.9 

7 

1 

1200 

190.7 

7 

1 

1800 

190.2 

7 

2 

000 

191.2 

7 

2 

600 

190.5 

7 

2 

1200 

190.6 

7 

2 

1800 

189.8 

7 

3 

000 

191.7 

7 

3 

600 

190.1 

7 

3 

1200 

189.9 

7 

3 

1800 

189.4 

7 

4 

000 

190.7 

7 

4 

600 

190.5 

7 

4 

1200 

193.5 

7 

4 

1800 

191.4 

7 

5 

000 

195.8 

7 

5 

600 

194.6 

7 

5 

1200 

193.1 

7 

5 

1800 

196.6 

7 

6 

000 

198.2 

7 

6 

600 

199.1 

7 

6 

1200 

196.3 

7 

6 

1800 

200.1 

7 

7 

000 

200.7 

7 

7 

600 

203.4 

7 

7 

1200 

205.4 

7 

7 

1800 

200.6 

7 

8 

000 

199.4 

7 

8 

600 

195.1 

sions  propagate  through  in  a  manner  similar  to  a 
chain-reaction  auto  collision.  A  possible  example 
of  such  an  event  is  at  about  2000  GMT  on  7  July 
(Julian  day  188)  in  lines  do  and  d,.  Note  the  sud¬ 
den  reversal  of  distance  in  do  at  this  time  after  con¬ 
vergence  had  been  occurring  for  some  time  in  line 
dj.  In  addition  to  this  simple  ridging  front,  simula¬ 
tions  have  shown  that  more  complex  propagating 
wave  effects  can  arise  if  the  coupling  between  the 
ice  strength  and  compactness  is  nonlinear. 

An  interesting  feature  in  the  recordings  was  the 
presence  of  occasional  sudden  slips  on  the  order  of 
tens  of  meters.  These  effects  probably  arise  due  to 
an  overshooting  phenomenon  in  the  ice  stress  that 
causes  the  distance  to  change  rapidly  and  then  re¬ 
main  approximately  at  a  constant  level.  An  exam¬ 
ple  is  shown  in  Figure  8,  where  a  time  series  of  the 
distance  do  has  been  drawn.  In  days  181  and  182, 
four  such  slips  occurred. 

A  perhaps  more  illustrative  example  of  fluctua¬ 
tions  is  shown  in  Figure  9,  which  shows  the  rela¬ 
tive  velocities  between  the  ship  (moored  to  a  floe) 
and  stations  O  and  J.  These  have  been  obtained 
first  through  central  differencing  the  raw  data  and 
then  by  filtering  with  a  low-pass  symmetric  filter 
with  241  weights  and  a  transition  band  from  48-  to 
60-min  periods  (Hibler,  1972).  The  relative  veloci¬ 
ties  are  very  large,  with  fluctuations  up  to  8  cm/s 
occurring.  Since  variations  in  the  ice  motion  will 
also  affect  current  water  strings  attached  to  ice 
floes,  these  results  emphasize  the  necessity  of 
making  detailed  ice  motion  measurements  to 
properly  interpret  current  measurements  made 
from  ice. 

The  stretch  speed  (i.e.  rate  of  change  of  dis¬ 
tance)  between  Polarbjdrn  and  the  on-ice  sites  was 
obtained  very  accurately  with  the  Del  Norte  sys¬ 
tem.  Figure  10  shows  the  spectra  of  the  stretch 
speeds  up  to  the  frequency  of  24  cycles  per  day. 
The  measurement  noise  lies  below  the  lowest  indi¬ 
cated  spectral  density  value.  For  sites  O  and  J, 
there  is  a  peak  close  to  2  cycles  per  day.  In  the 
higher  frequencies  the  spectral  density  falls  off 
with  a  red  noise  character.  A  comparison  with  the 
radar  observations  (Fig.  11)  shows  that  the  radar 
also  gives  the  peak  at  about  2  cycles  per  day,  but 
then  the  spectral  density  falls  off  much  more 
rapidly  than  in  the  case  of  Del  Norte  data.  Due  to 
its  inaccuracy  it  is  unlikely  that  the  radar  can  re¬ 
solve  higher  frequencies  than  about  3  cycles  per 
day. 

Overall  this  data  set  represents  a  unique  set  of 
measurements  of  the  deformation  field  of  a  series 
of  interacting  floes  floating  in  the  ocean.  Because 
of  the  fine  temporal  resolution,  accuracy  of  the 
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Table  3.  Linear  first  order  fit  for  the  Del  Norte  data  (6-hr  displacements)  in  north-east  coordinate 
system. 


Oata 

Mo 

Day 

Hr 

Oeu 

Dxu 

Oyu 

Dov 

Dxv 

Dyv 

El 

Ell 

Ei 

(2 

*i 

Vo r 

Fit 

«*T) 

(cm/%) 

(l/H) 

(I/M) 

(ca/s) 

(l/H) 

(l/H) 

(l/H) 

(l/H) 

(l/H) 

(l/H) 

(dag) 

(0/H) 

(i) 

6 

29 

0 

0.4  1 

1.90 

0.50 

1.38 

2.64 

-2.20 

-0.30 

5.16 

2.43 

-2.73 

57.2 

0.61 

91.8 

6 

29 

600 

•1.31 

-0.49 

-1.19 

-2.13 

-1.50 

-0.93 

-1.42 

2.81 

0.69 

-2.11 

-26.6 

-0.11 

64.7 

6 

29 

1200 

-0.01 

-0.46 

-0.07 

0.01 

0.37 

0.14 

-0.32 

0.66 

0.17 

-0.49 

40.0 

0.13 

99.9 

6 

29 

1000 

-0.08 

0.53 

-0.76 

-0.12 

0.34 

-0.38 

0.14 

1.00 

0.57 

-0.43 

-69,9 

0.32 

98.0 

6 

JO 

0 

-0.21 

0.44 

-0.30 

-0.40 

-1.40 

-0.62 

-0.17 

2.00 

0.91 

-1.09 

-47.1 

-0.32 

93.7 

6 

30 

600 

-0.01 

-0.11 

0.06 

-0.04 

-0.24 

0.03 

-0.07 

0.22 

0.08 

-0.15 

-40.8 

-0.09 

98.5 

6 

30 

1200 

0.21 

-0.53 

0.65 

0.24 

0.25 

0.23 

-0.30 

1.18 

0.44 

-0.74 

44.5 

-0.12 

91.4 

6 

30 

1000 

-0.13 

0.25 

-0.60 

-0.11 

-0.66 

-0.02 

0.22 

1.29 

0.76 

-0.53 

-50.2 

-0.02 

97.0 

7 

1 

0 

0.06 

-0.13 

0.36 

0.09 

0.24 

-0.04 

-0.17 

0.61 

0.22 

-0.39 

35.9 

-0.04 

96.6 

7 

1 

600 

-0.10 

-0.29 

0.32 

-0.06 

-0.42 

0.12 

-0.17 

0.43 

0.13 

-0.30 

-69.1 

-0.21 
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Fig.  9.  Rate  of  change  of  the  distances  d>  and  d,  during  30 
hours  in  late  part  of  the  experiment.  The  time  series  has 
been  filtered  with  a  low-pass  filter  ( transition  band  48-  to 
60-minute  periods). 
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Fig.  11.  The  spectra  of  the  rate  of  change  of .  he  distance  d,  as  measured 
with  the  Del  Norte  system  and  by  using  the  ship’s  radar. 


measurements,  and  coincidence  of  current  and 
wind  measurements,  this  data  should  greatly  aid  in 
sorting  out  the  physics  of  marginal  ice  zone  ice  dy¬ 
namics  and  kinematics.  Such  efforts  are  currently 
in  progress. 
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Stresses  due  to  ice/ice  interactions  have  been  modeled  to  obey  a  plastic  constitu¬ 
tive  law  (Hibler,  1979)  or  a  Newtonian  viscous  constitutive  law  (Kulakov  et  al., 
1979).  The  applications  of  these  models  have  successfully  explained  long-term,  over¬ 
all  ice  motion  in  the  Greenland  Sea  and  the  Baltic  Sea  (e.g.  Hibler,  1979;  Lepparan- 
ta,  1981). 

However,  near  the  ice  edge  floe  characteristics  are  different  from  the  interior  pack 
ice.  As  a  consequence,  while  observed  mesoscale  ice  edge  phenomena  [for  instance, 
the  occurrence  of  ice-edge  jets  (Johannessen  et  al.,  1983)]  can  be  reproduced  by  ideal 
plastic  constitutive  laws  through  transient  adjustment  states  (Lepparanta  and  Hib¬ 
ler,  1985),  it  is  desirable  to  understand  ice  rheology  in  the  marginal  ice  zone  from  a 
more  fundamental  point  of  view. 

During  MIZEX  83,  a  considerable  amount  of  precision  ice  deformation  data 
(Hibler  and  Lepparanta,  1984)  were  obtained.  These  data  provide  a  basis  for  sorting 
out  the  physics  that  governs  the  ice  rheology  in  the  MIZ. 

The  deformation  data,  together  with  other  MIZEX  83  observations,  show  that  in 
the  MIZ  the  moving  ice  field  consists  of  a  granular  material  with  various  sizes, 
shapes,  and  thicknesses.  This  material  undergoes  a  deformation  that  is  quite  chaotic 
when  a  fine  temporal-resolution  record  is  studied.  Ice  floes  fluctuate  from  their 
mean  positions  and  collide  with  neighboring  ones  without  forming  ridges.  The  fre¬ 
quency  of  these  occurrences  is  in  the  order  of  10/day.  Assuming  the  average  floe  in 
the  MIZ  has  a  diameter  of  50  m  and  a  thickness  of  1  m,  and  the  ice  compactness  is 
about  0.8,  the  collisional  momentum  transfer  rate  is  then  of  the  order  of  1  dyne/cm. 
This  magnitude  is  comparable  with  other  forcing  mechanisms  that  exist  in  the  field. 
Hence,  it  is  of  considerable  interest  to  analyze  the  internal  stresses  as  a  result  of  the 
collisional  ice  interactions. 

In  this  note  we  outline  a  preliminary  theory  for  determining  the  rheology  of  a 
broken  icefield  due  to  floe  collisions.  The  theory,  which  is  based  on  an  idealized 
granular  flow  concept,  quantifies  the  deformation-induced  velocity  fluctuations  as 
well  as  the  internal  stresses.  Initial  comparison  of  the  theoretical  velocity  fluctua¬ 
tions  are  then  made  with  the  Held  data  obtained  in  MIZEX  83. 


Consider  two  disks,  A  and  B,  embedded  in  a  deforming  ideal  granular  material,  as 
shown  in  Figure  1 .  When  these  two  disks  collide,  the  amount  of  momentum  transfer 
depends  on  their  relative  velocity  and  their  material  properties  as  characterized  by 
density  a,  the  restitution  coefficient  t,  and  the  frictional  coefficient  /<. 

For  a  given  deformation  tensor 


du/dx  du/dy 


dv/dx  dv/dy 


the  relative  velocity  between  A  and  B  is 


Va"  =  d(4t  cos9+  ~iysine)i+ 

D  (~y  cosd  +  -^ysin^jj  +  \fl  v'  (2) 

where  v'  is  an  isotropic  random  fluctuation  coi-.ponent  of  a  single  disk.  After  the 
collision,  from  linear  momentum  conservation  and  energy  loss  due  to  friction  and 
incomplete  restitution,  the  momentum  transfer  can  be  obtained  as 

AMg  cos30-/jcos20  sin0)^  +  (sin0  cos2 9 -fJ sin2 6  cos0)^  + 

+  (cos0  sin20-psin30)|^j/  +  j(sin0  cos20  +  pcos30) 

+  (sin20cos0+psin0cos20)  +  (sin30  +  #icos0  sinJ0)|^j/J 


+,ffi  \/2  v  sina N 


After  statistically  averaging  all  possible  collision  arrangements  taking  place  between 
A  and  B, 

for  case  I  in  Figure  1 ,  and 

<s> 

for  case  2  in  Figure  1 . 

The  stress  tensor  generated  by  these  collisions  is  modeled  as  the  rate  of  total  mo¬ 
mentum  transfer  across  a  surface  defined  by  all  disks  B  that  occupy  the  unit  control 
surface.  Hence 

Tii  =  P.-/-MT,  (6) 


D 


)  O  (B. 

go. 


h/*\? 


Fig.  1.  A  deforming  field  of  disks  and  two  blow-up  cases  for  a  given  deformation 
tensor. 


where  p ,  =  number  of  disks  on  a  unit  surface  =  p  (if  isotropic  distribution) 

/  =  v  '/2S;  S  =  average  gap  size  (7) 

AMj  =  Xj  component  of  Mf, 

Substituting  (4)  and  (5)  into  (6),  we  obtain 

r  =  ms~  -f-p  {(1+2pX7>{>)/*(1-m)Z>]-^v'/|  (8) 

The  fluctuation  component  v '  has  to  do  with  the  energetics  of  the  deforming  granu¬ 
lar  material.  Stronger  deformation  induces  more  frequent  collisions,  hence  a  higher 
v'.  Higher  v'  increases  the  frequency  of  collisions.  Each  collision  involves  three 
basic  dissipation  mechanisms:  that  due  to  moving  through  the  surrounding  fluid, 
that  due  to  an  incomplete  restitution  after  a  collision,  and  that  due  to  the  frictional 
force  between  disks  during  contact.  These  dissipations  damp  out  collisional  energy 
and  reduce  v Assuming  the  time  needed  to  adjust  the  value  of  v '  to  a  gradually 
changing  deformation  field  is  short,  balance  between  the  deformational  increase  and 
the  collisional  damping  of  v '  can  be  arrived  at  instantaneously.  Therefore, 

fP  =  ^  •/•£  (9) 

where  Q  is  the  strain  tensor  and  E  is  the  average  energy  loss  per  collision.  The  energy 
loss  due  to  a  collision  shown  in  Figure  1  is 

•/2  ms  [  VA*  +  -  Vf1-  Vr\  +SCD  Q/Ddv  '*  (10) 

where  *  denotes  post-collision  conditions.  Again  averaging  throughout  all  possible 
collision  arrangements,  we  obtain 

*■  ¥  [(.-«=)-»“<'*«, •!  NS*f )'  *  '*1-1  f/ 

+  w  +  +  1  £p?  El\  (11) 

*  D  l  4  ir  4  it  D  psl 

Substituting  (8)  and  (11)  into  (9),  we  obtain  an  equation  for  v '  as  the  following: 

k,v ,}  +  k,v '  +  k,  =  0  (12) 


» 


t  -Jl-e*  Ml-t-e)  ifjl+e)*  4  CDS  pf 
1  4  it  4  it  D  Ps 


k3  =^fl-eJ-M10+e)1](J/,/?-,/l/2)-(l+e)37  [(2+^)/f - 2(l-#i)/} ] 

.  _  bu  bv  bu  bv  ,,/9m  3v\J 

Il=  bx  by ’’  /j  =  Ax  3>-  "  /4\d>-  9x/ 


(13) 


The  above  quadratic  equation  for  v '  can  be  solved  if  the  material  properties  Q/>  Q»  «. 
M,  and  the  fluid  drag  coefficient  CD  are  given,  and  the  disk  size  D,  average  gap  size  S, 
and  deformation  invariants  /,  and  h  are  provided. 

For  MIZEX  83  ice  data,  the  material  properties  are  estimated  as  q,/q,  =  0.9,  e  = 
0.2,  n  =  0.2,  CD  =  0.5,  D  =  50  m,  and  S  =  5  m  (which  corresponds  to  a  compact¬ 
ness  of  0.8).  The  expected  v '  from  the  theory  is 

v.'=  DI-0.39/,  +  V  0.69/M).  25/i  |.  (14) 


Since  /  =  v  '/2S,  the  stress  tensor  may  be  obtained  when  the  above  is  substituted 
into  (8). 

During  the  MIZEX  83  mesoscale  ice  dynamics  program,  the  positions  of  four 
transponders  relative  to  the  ship  were  recorded  at  3-min  intervals  (Hibler  and  Leppa- 
ranta,  1984).  The  data  were  then  filtered  through  a  low-pass  filter  of  transition  band 
48-60  minutes.  Adopting  the  Del  Norte  coordinate  system,  the  velocity  field  was 
then  computed  and  resampled  at  half-hour  frequency.  This  velocity  field  is  denoted 
here  as  V.).  Because  the  ship  was  used  as  a  reference  point  in  the  Del  Norte 
system  V.  of  the  ship  is  zero. 

A  linear  first-order  fit  procedure  is  then  adopted  to  obtain  the  ship’s  velocity  ( K) 
and  the  components  in  the  deformation  tensor  (du/dx,  du/d y,  d  v/dx,  d  v/d y).  This 
deformation  tensor  and  the  filtered  position  data  for  the  four  transponders  are  then 
utilized  to  give  the  strain-predicted  velocity  at  the  transponders’  and  the  ship’s  loca¬ 
tions.  These  velocity  data  are  denoted  here  as  /„(Ff/).  These  strain-predicted 
velocities  are  viewed  as  the  mean  velocities  of  the  four  transponders. 

Let 


y.  =  u.i  +  vj 
V,  =  u.i+vj 

v„  =  u„+  V  (15) 

in  the  Lagrangian  coordinate  system.  The  actual  velocities  of  the  transponders  and 
the  ship  are 

V}  =  V.-V,  =  {u.-u,)i  +  (v.-v.y  for  the  transponders 


Vt  =  -K  •  -uj-vj 


for  the  ship 


(16) 
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Fig.  2.  Comparison  of  fluctuations. 
v, '  from  file  EXNFIL,  v„ '  from  file 
MEANMY.  A  conversion  factor  of 
0.014  (cm/sec)/ [cm(%/hr)J  has  to 
be  multiplied  to  v,  ’  to  put  it  into  the 
same  unit  ui  v, 


and  the  linear  fit  or  the  mean  velocities  of  the  transponders  and  the  ship  are 
V\f  =  V,f-V,  =  (ulf-u,)i  +  (vt/-v,)j  for  the  transponders 

Vlf  =  0  for  the  ship.  (17) 

The  velocity  fluctuation  from  measured  data  is  obtained  as 

v:  =  V,  L  \Vt-V)f\.  (18) 

The  above  summation  is  taken  over  the  four  transponders  and  the  ship. 

The  two  velocity  fluctuation  data  sets  v.'and  vl  are  computed  from  (14)  and  (18) 
at  half-hour  intervals.  These  two  data  sets  are  smoothed  with  an  8-hr  block  average 
procedure.  The  data  are  then  represented  in  histograms  at  an  8-hr  step  size  and  the 
correlation  coefficient  between  the  two  histograms  is  obtained.  A  typical  histogram 
presentation  of  the  correlation  between  v.'and  is  given  in  Figure  2. 

Since  the  time  series  of  the  data  have  wide  spectra,  correlation  between  v.'and  v^  is 
investigated  under  six  different  cases.  Let 

,>.0(data)  (19) 

represent  the  low-pass  filtered  data  with  band  width  from  P  to  Q  hr.  Three  filtered 
data  sets  Fj  j-iof/i./i).  ^.$-io(KJ),  and  F2.j_io(  FJy)  are  constructed.  Five  velocity  fluctu¬ 
ation  files  are  obtained  with  or  without  the  filtering  procedure.  Table  1  illustrates 
how  these  five  files  are  obtained. 

The  result  of  the  correlation  and  the  mean  and  standard  deviation  of  pairs  of  his¬ 
tograms  for  v.'and  vj,  are  listed  in  Table  2,  where  D  =  50  m,  q./q,  =  0.9,  t  =  a  = 
0.2,  and  Co  =  0.5  are  used  in  (14). 
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Table  1.  Files  of  velocity  fluctuation. 


» 


File 

name 

Data 

provided 

Data  used 

EXNF1L 

(eq  14) 

/,  no  filter 

EXYFIL 

(eq  14) 

Fis-vVM 

MEANMN 

(eq  18) 

K‘.  y\,  no  filter 

MEANMY 

(eq  IS) 

Kao  filter,  F2.s-io(Kf/) 

MEAYMY 

(eq  18) 

*2.5-IO(*?.  »?y» 

Table  2.  Correlation  of  velocity  fluctuations. 


Theory 

Mean 

S.D. 

Field 

data 

Mean 

S.D 

Correlation 

coefficient 

EXNF1L 

0.01 

0.007 

MEANMN 

0.60 

0.40 

0.9 

MEANMY 

1.19 

0.77 

0.83 

MEAYMY 

0.35 

0.28 

0.9 

EXYFIL 

0.008 

0.006 

MEANMN 

0.72 

MEANMY 

0.59 

MEAYMY 

0.89 
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The  purpose  of  this  note  is  to  demonstrate  the  overall  procedure  for  determining 
the  rheology  of  a  broken  ice  field  due  to  floe  collisions.  However,  for  a  consistent 
theory  some  modifications  are  needed.  Most  notably  these  equations  have  not  in¬ 
cluded  the  mean  flow  in  the  determination  of  the  rheology,  which  is  particularly  im¬ 
portant  in  divergence.  Also,  although  good  correlation  is  obtained  with  the  observed 
deformation  data  the  observed  fluctuations  are  larger  than  predicted.  A  more  com¬ 
plete  analysis  addressing  these  issues  is  currently  underway  and  will  be  submitted  for 
journal  publication  in  the  near  future. 
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INTRODUCTION 

The  summer  decay  and  retreat  of  the  marginal 
ice  zone  is  controlled  by  a  variety  of  complex  in¬ 
teractions  taking  place  among  the  ice,  ocean  and 
atmosphere.  Melting  alters  not  only  the  properties 
of  the  ice  cover,  but  also  the  structure  of  the  upper 
ocean.  Horizontal  gradients  in  meltwater  input 
and  the  proximity  of  warmer  and  more  saline 
water  can  give  rise  to  a  variety  of  mesoscale  phe¬ 
nomena  in  the  ocean,  while  floe  breakup  and  in¬ 
creasing  stratification  beneath  the  ice  affect  the  re¬ 
sponse  of  the  ice  to  winds  and  currents  and  the 
rate  at  which  heat  is  transferred  from  the  water  to 
the  ice.  Both  the  atmosphere  and  ocean  supply 
energy  needed  to  drive  the  melt.  The  primary 
source  of  energy  from  the  atmosphere  is  short¬ 
wave  radiation,  and  the  properties  of  the  ice  which 
affect  its  absorption  and  distribution  must  be  con¬ 
sidered  as  important  variables  in  the  system.  Heat 
transfer  from  the  water  to  the  ice  occurs  at  both 
the  edges  and  bottoms  of  floes  and  is  a  major  fac¬ 
tor  in  the  summer  decay  cycle.  Much  of  the  heat 
supplied  to  the  ice  by  the  water  appears  to  be  ab¬ 
sorbed  shortwave  radiation  which  enters  the  upper 
ocean  through  leads  and  areas  of  thin  ice.  How¬ 
ever,  the  rate  at  which  this  heat  is  lost  to  the  ice  is 
complicated  by  lateral  advection  of  ice  and  water 
across  the  edge,  mixed  layer  stratification,  waves, 
floe  interactions  and  mesoscale  circulations  asso¬ 
ciated  with  the  MIZ.  Quantitative  information  on 
the  role  of  these  processes  is  meager. 

At  present  we  are  unable  to  answer  many  basic 
questions  regarding  the  heat  and  mass  balance  of 
the  summer  MIZ.  For  example,  what  is  the  relative 
importance  of  the  atmosphere  and  ocean  in  the  de¬ 
cay  cycle?  What  fraction  of  the  mass  loss  occurs  at 
the  extreme  ice  edge?  How  is  solar  energy  ab¬ 


sorbed  in  the  water  partitioned  between  lateral 
melting,  bottom  ablation,  temporary  storage  in 
the  ocean,  and  loss  to  the  atmosphere?  What  are 
the  dominant  processes  involved  in  the  transport 
of  heat  from  the  water  to  the  ice?  To  what  extent 
does  freshwater  input  from  the  melting  ice  retard 
the  vertical  exchange  of  heat  with  the  deeper 
ocean?  How  does  ice  movement  and  proximity  to 
the  edge  affect  heat  transport  and  density  struc¬ 
ture  in  the  summer  mixed  layer?  To  answer  such 
questions  requires  concurrent  data  on  the  state  of 
the  ice  (concentration,  thickness,  floe  size  distri¬ 
bution,  snow  and  pond  coverage),  mass  changes  at 
floe  boundaries,  energy  input  to  the  ice  and  ocean, 
ice  movement,  ocean  structure,  and  heat  exchange 
processes  between  system  components.  Although 
the  program  outlined  below  is  concerned  primarily 
with  temporal  changes  in  the  heat  and  mass  bal¬ 
ance  of  the  ice  pack,  it  should  provide  information 
needed  for  at  least  a  preliminary  examination  of 
many  of  the  above  questions.  Because  shortwave 
radiation  is  a  principal  agent  in  the  decay  cycle, 
the  work  has  a  strong  focus  on  properties  and 
processes  that  affect  its  interaction  with  the  ice 
and  upper  ocean. 

EXPERIMENTAL  PROGRAM 

Our  overall  plan  for  MIZEX  83  was  to  instru¬ 
ment  a  few  typical  floes  in  the  MIZ  and  to  monitor 
mass  changes  occurring  at  the  top,  bottom  and 
sides  of  these  floes.  Other  measurements  included 
incoming  radiation,  ice  thickness,  snow  depth, 
pond  coverage,  internal  ice  temperatures,  radia¬ 
tive  properties,  and  heat  content  changes  in  leads. 


We  plan  to  use  these  data  to  determine  the  relative 
importance  of  the  various  energy  sources  and  to 
construct  local  heat  and  mass  budgets  for  the  indi¬ 
vidual  floes.  With  these  local  budgets  and  esti¬ 
mates  of  ice  concentration  and  pond  coverage  de¬ 
rived  from  aerial  photography  it  should  ultimately 
be  possible  to  obtain  area-integrated  totals  for  the 
region.  Of  particular  concern  in  the  study  was  the 
input  of  solar  energy  to  the  ocean  and  the  factors 
that  control  its  transfer  to  the  ice.  Analysis  of  the 
above  data  in  conjunction  with  results  from  the 
oceanographic  program  should  yield  detailed  in¬ 
formation  on  the  role  of  the  upper  ocean  in  the  de¬ 
cay  and  retreat  of  the  ice  and  allow  us  to  gauge  its 
importance  in  relation  to  atmosphere-ice  interac¬ 
tions. 

Heat  and  mass  balance  observations  during  the 
drift  phase  were  concentrated  on  two  floes 
(“Tent”  Floe  and  “Spot”  Floe)  located  about  150 
m  south  of  the  Polar bj dr n.  Tent  Floe,  measuring 
about  25  x75  m  in  the  horizontal,  averaged  1.9  m 
in  thickness  and  was  composed  largely  of  first- 
year  ice.  The  somewhat  thicker  (2.4  m)  Spot  Floe 
measured  about  50  x75  m  and  appeared  to  contain 
areas  of  both  first-year  and  second-year  ice.  Elev¬ 
en  hot  wire  thickness  gauges  were  installed  in  these 
two  floes  to  monitor  mass  changes  at  the  bottom 
of  the  ice.  An  additional  four  gauges  were  placed 
in  the  ice  next  to  the  ship  to  see  to  what  extent 
waste  heat  from  the  ship  might  affect  the  melting. 
Surface  ablation  was  measured  at  27  sites  on  these 
three  floes.  Wall  profiles  and  observations  of  edge 
erosion  were  made  at  several  locations  around 
Tent  Floe.  Salinity  and  temperature  profiles, 
taken  in  leads  during  the  drift  and  ice  edge  phases, 
were  used  to  calculate  the  heat  content  of  the  up¬ 
per  several  meters  of  the  water  column.  Total  in¬ 
coming  shortwave  and  longwave  radiation  fluxes 
were  monitored  continuously  over  the  ice  on  Tent 
Floe  throughout  the  drift  period.  Complementary 
data  on  the  effect  of  cloudiness  on  the  spectral  dis¬ 
tribution  of  incident  irradiance  were  also  ob¬ 
tained.  Spectral  and  total  albedo  measurements 
were  made  over  bare  multiyear  ice,  snow,  ponds, 
and  thin  first-year  ice.  Conductive  heat  fluxes  in 
the  ice  were  estimated  from  internal  temperature 
data  in  Spot  Floe. 


RESULTS 

Conditions  during  the  drift  phase  were  unex¬ 
pectedly  quiescent— ice  concentration  was  high, 
leads  and  floes  were  small,  wave  action  was  mini¬ 
mal  and  there  was  little  in  the  way  of  floe  breakup 


or  mechanical  erosion.  Because  the  drift  took 
place  early  in  the  melt  season,  albedos  were  still 
high,  slowing  the  surface  melt  and  producing  only 
a  slight  freshening  of  the  upper  ocean.  Ice  decay 
was  probably  close  to  the  minimum  amount  that 
can  be  expected  in  the  summer  MIZ.  Many  floes  in 
the  region  contained  significant  amounts  of  algae 
and  dirt,  lowering  their  albedos  and  increasing  the 
amount  of  shortwave  energy  absorbed.  Offsetting 
this  effect,  however,  was  the  small  melt  pond  cov¬ 
erage  on  the  thicker  (>  1  m)  floes,  which  caused 
them  to  have  somewhat  higher  albedos  than  simi¬ 
lar  ice  in  the  central  pack  or  in  coastal  areas.  Floes 
thinner  than  1  m  tended  to  be  quite  smooth  and 
were  frequently  covered  almost  entirely  by  shallow 
ponds.  Some  large  multiyear  floes  with  pond  cov¬ 
erage  typical  of  the  central  Arctic  were  sighted 
during  aircraft  overflights,  but  were  unusual.  The 
great  majority  of  floes  in  the  vicinity  of  the  ship 
appeared  to  be  similar  in  character  to  Spot  and 
Tent  Floes. 

Bottom  ablation 

Daily  ablation  totals  at  each  of  the  thickness 
gauge  sites  are  given  in  Table  1.  Cumulative  bot¬ 
tom  ablation  over  the  experimental  period  is 
sketched  in  Figure  1.  Average  melt  rates  varied  be¬ 
tween  5  and  9  mm/day.  Standard  deviations  were 
large.  Part  of  the  explanation  for  the  large  devia¬ 
tions  was  a  definite  tendency  for  thicker  ice  to  ab¬ 
late  more  rapidly.  Figure  2  shows  the  average  melt 
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Fig.  I.  Cumulative  mass  losses  oc¬ 
curring  at  the  top  (solid  curve )  and 
bottom  (dashed  curve )  of  floes  near 
the  Polarbjdrn  during  the  drift 
phase. 


Tabic  1.  Bottom  ablation  (mm/day)  on  floes  near  the  Polarbjprn  during  the  drift  phase.  Gauges  D1-D3  were  lo¬ 
cated  on  “Dirty  Floe’’  adjacent  to  the  ship;  gauges  T1-T6  were  located  on  “Tent  Floe,"  the  site  of  the  mixed  layer  studies; 
and  gauges  S1-S5  were  located  on  “Spot  Floe.’’  h,  was  the  initial  ice  thickness  when  the  gauges  were  installed.  The  nega¬ 
tive  value  at  S5  on  5  July  is  the  result  of  meltwater  input  beneath  the  ice  causing  formation  of  a  false  bottom. 
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Fig.  2.  Average  melt  rate  at  the  bottom  of 
the  ice  vs  initial  ice  thickness  during  the 
period  3-8  July  1983. 


rate  (f)  over  the  final  six  days  of  the  drift  plotted 
against  the  initial  ice  thickness  (h0).  The  solid  line 
in  Figure  2  represents  a  least  squares  fit  to  the 
data:  /  =  -7.64  +  6.9/j0  ( R 2  =  0.83).  A  strong  de¬ 
pendence  of  /  on  h  was  not  expected  over  such  a 
small  thickness  range  and  we  are  unsure  as  to  the 
reason.  The  two  largest  values  of/ were  measured 
in  small  ridge  keels,  but  thickness  measurements 
were  not  sufficiently  dense  to  provide  a  good  defi¬ 
nition  of  bottom  roughness.  Despite  the  relatively 
small  amount  of  bottom  melting,  ablation  rates 
were  still  much  larger  than  those  observed  in  areas 
of  perennial  ice. 

Snow  cover  and  surface  ablation 

Although  melting  at  the  upper  surface  was  al¬ 
ready  in  progress,  a  substantial  amount  of  snow 
still  covered  the  ice  at  the  beginning  of  the  drift. 
Measurements  from  a  121-point  survey  grid  (Ta¬ 
ble  2)  showed  snow  thicknesses  varying  between  0 
and  123  cm,  with  an  average  of  43.3  cm.  Snow 
density  remained  essentially  constant  throughout 
the  period  at  400  kg/m’.  Surface  ablation  totaled 
about  115  kg/m!  during  the  period  (Figure  1),  a 
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Fig.  3.  Histogram  of  snow  depths  near  Polar- 
bj0rn,  27  June  1983. 

mass  loss  roughly  double  that  observed  on  the  un¬ 
derside  of  the  ice.  Daily  ablation  amounts  are  tab¬ 
ulated  for  each  of  these  sites  in  Table  3.  Area- 
weighted  surface  melt  rates  for  bare  ice,  snow- 
covered  ice,  and  the  combined  average  are  present¬ 
ed  in  Figure  4.  The  decreasing  rate  of  meltwater 
generation  by  the  snow  and  the  corresponding  in¬ 
crease  by  the  ice  simply  reflect  changes  in  the  rela¬ 
tive  area  covered  by  snow  and  bare  ice. 

The  heavy  snow  cover  encountered  in  the  MIZ 
was  in  striking  contrast  to  regions  of  perennial  ice 
and  to  coastal  portions  of  the  Arctic  Basin  where 
snow  typically  vanishes  by  1  July.  Extensive  snow 
cover  was  also  noted  three  weeks  later  in  the 
Greenland  Sea  during  the  transect  of  the  Polar- 
stern  to  the  Greenland  coast.  At  this  point  it  is  dif¬ 
ficult  to  be  certain  whether  lower  energy  input  at 
the  surface  or  greater  snowfall  was  primarily  re¬ 
sponsible  for  the  large  amounts  of  snow,  but  we 


Fig.  4.  Area-weighted  melt  rate  averages  at 
the  surface  of  the  ice  cover  during  the  drift 
phase.  Solid  curve  represents  total  rate  of  mass 
loss,  while  the  dashed  and  dotted  curves  represent 
the  respective  contributions  from  ice  and  snow  to 
the  total. 

suspect  the  latter.  The  significance  of  the  extensive 
snow  cover  is  that  it  reduces  the  absorption  of 
solar  radiation  and  protects  the  surface  of  the  ice 
from  erosion.  This  slows  the  decay  cycle  and 
means  that  heat  from  the  water  must  be  a  major 
factor  in  the  retreat  of  the  MIZ  during  early  sum¬ 
mer. 


Table  2.  Snow  depth  (mm)  survey  taken  on  rough  2nd  year  floe  near  Polar- 
bjdrn,  27  June  1983.  Grid  spacing  was  3  m.  Melting  was  in  progress  and  the  snow  was 
nearly  isothermal.  Mean  snow  depth  was  433  mm  with  a  standard  deviation  of  274  mm. 
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Table  3.  Summary  of  surface  ablation  (mm)  measurements  made  during  the  drift  phase.  Daily  mean  mass  loss¬ 
es  are  given  in  kg/m*  and  specific  loss  rates  In  kg/m'  day.  A  “/”  beside  a  number  indicates  a  transition  from  snow- 
covered  to  bare  ice,  values  to  the  right  of  the  slash  representing  ice  ablation  and  those  to  the  left  snow  ablation. 
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Fig.  5.  Ice  wall  profiles  taken  on  the  edge  of  Tent  Floe  on  5 
and  7  July. 


Lateral  erosion 

Ice  wall  profiles  measured  on  Tent  Floe  were  of 
two  general  types:  1)  those  shaped  primarily  by 
thermal  processes,  which  had  an  underwater 
ledge,  frequently  accompanied  by  an  overhanging 
lip  above  the  waterline,  and  2)  those  shaped  by 
mechanical  erosion  from  floe  interactions,  which 
were  more  vertical  and  lacked  the  underwater 
bulge.  Undercutting  on  the  MIZEX  floes  was  sub¬ 
stantially  less  than  the  amount  we  observed  during 
July  1982  in  the  dynamically  active  MIZ  west  of 
Prince  Patrick  Island.  Figure  5  shows  two  profiles 
taken  on  5  and  7  July.  Despite  the  small  wave  ac¬ 
tivity,  the  rate  of  mass  loss  was  large  when  com¬ 
pared  to  those  at  the  top  and  bottom  of  the  floe. 
Horizontal  retreat  of  the  vertical  ice  wall  below 
the  water  was  about  S  cm/day,  while  retreat  of  the 
sloping  underwater  “ramp”  averaged  20-30  cm/ 
day.  Horizontal  retreat  of  the  notch  at  the  water¬ 
line  was  12  cm/day,  roughly  the  same  as  that  of 
the  overhanging  lip.  At  the  start  of  the  observa¬ 
tional  period,  the  leads  contained  a  great  deal  of 
ice  debris,  which  was  probably  generated  as  a  re¬ 
sult  of  floe  jostling  by  the  ship.  This  debris  melted 
rapidly  during  the  first  two  days.  Little  debris  was 
sighted  after  this  period,  indicating  that  mechani¬ 
cal  erosion  on  floe  edges  was  minimal. 

Melt  pond  and  bare  ice  coverage 

Melt  pond  coverage  was  much  less  than  that  ob¬ 
served  in  the  central  Arctic  or  in  coastal  regions. 
Pond  area  reached  a  maximum  of  about  20 %  for 
1-2  days,  then  decreased  rapidly  to  about  5%  by  8 
July  (Figure  6).  The  reason  for  the  smaller  cover¬ 
age  was  the  smaller  size  of  the  floes,  which  made  it 
relatively  easy  for  meltwater  to  find  a  path  to  the 


ocean.  An  exception  to  this  was  the  thin  (<  1  m) 
floes,  which  were  flat,  and  whose  surfaces  were 
close  to  the  waterline.  Such  floes  tended  to  be  cov¬ 
ered  by  extensive  but  very  shallow  ponds.  Since 
the  more  heavily  ponded  multiyear  floes  were  not 
a  significant  factor,  surface  meltwater  storage  was 
low,  and  freshwater  input  to  the  ocean  in  the  drift 
area  should  be  closely  correlated  with  measured 
ablation. 

Figure  6  also  shows  the  relative  area  covered  by 
bare  ice,  based  on  data  from  the  surface  ablation 
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Fig.  6.  Relative  area  covered  by 
ponds  and  bare  ice  during  the  drift 
phase. 
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array.  Because  stakes  were  frequently  located  in 
higher  areas,  these  values  may  be  somewhat  biased 
toward  areas  of  bare  ice.  It  should  be  remem¬ 
bered,  however,  that  values  for  pond  and  bare  ice 
areas  apply  only  to  the  three  instrumented  floes. 
Aerial  photos  in  the  vicinity  of  the  ship  disclose 
numerous  thin  pond-covered  floes,  so  when  re¬ 
gional  totals  are  ultimately  obtained,  they  will 
probably  indicate  less  snow  and  greater  pond  cov¬ 
erage. 

Internal  ice  temperatures 

To  determine  if  there  was  any  significant  con¬ 


duction  of  heat  within  the  ice,  temperatures  were 
measured  at  five  levels  in  3-m  ice  near  a  thickness 
gauge  in  the  middle  of  Spot  Floe.  The  results 
(Table  4)  indicate  that  temperatures  in  the  lower 
meter  of  the  ice  were  essentially  isothermal.  Tem¬ 
peratures  in  the  upper  portion  of  the  ice  tended  to 
be  more  variable,  reflecting  differences  in  cloudi¬ 
ness  and  air  temperature  on  the  different  days.  Al¬ 
though  there  was  slight  warming  over  the  nine-day 
period,  heat  conduction  near  the  bottom  of  the  ice 
was  small,  and  bottom  ablation  rates  should  pro¬ 
vide  a  direct  measure  of  the  oceanic  heat  flux. 


Table  4.  Internal  Ice  temperatures  (°C)  measured  at  the  center 
of  Spot  Floe  during  the  drift  phase.  Initial  ke  thickness  was  2.97 


Fig.  7.  Comparison  of  salinity  profiles  taken  in  leads  at  three  sites  dur¬ 
ing  the  drift  and  ice  edge  phases.  Site  A  (dotted  curve)  was  a  25-m-wide 
lead  adjacent  to  Spot  Floe  on  3  July;  Site  B  was  a  3-  to  4-m-wide  lead  in  a  re¬ 
gion  of  30%  ice  concentration  on  21  July  (dashed  curve):  and  Site  C  (solid 
curve)  was  on  the  windward  side  of  a  50-m-wide  lead  in  a  region  of  low  ice 
concentration  on  22  July. 
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Heat  content  of  leads 

To  obtain  information  on  the  distribution  of 
heat  in  leads,  CTD  profiles  were  taken  in  the  upper 
2-4  m  of  the  water  column  using  hand-held  Sea¬ 
bird  temperature  and  conductivity  cells.  Results 
are  summarized  in  Figures  7-9.  Measurements 
during  the  drift  phase  were  made  on  3  July  in  a 
25-m-wide  rectangular  area  of  open  water  adja¬ 
cent  to  Spot  and  Tent  Floes.  Profiles  were  taken  in 
the  center  of  the  lead,  as  well  as  on  the  upwind  and 
downwind  sides.  Although  the  heat  content  of  the 
lead  was  about  15%  larger  on  the  downwind  side, 


profiles  were  similar  at  all  three  sites,  showing  a 
lower  salinity  meltwater  layer  in  the  upper  10-20 
cm.  There  was  no  evidence  of  meltwater  entrain¬ 
ment  below  this  level.  Ice  concentration  at  this 
time  was  about  85%,  and  the  cumulative  melting, 
averaged  over  the  region,  totaled  about  60  kg/m!. 
Heat  content  in  the  lead  was  20-25%  of  the  daily 
solar  input,  indicating  rapid  transport  of  heat  to 
the  walls  and/or  under  the  ice. 

CTD  data  were  also  gathered  in  several  leads 
further  south  (79  °30 ')  during  the  ice  edge  phase. 
On  21  July  six  profiles  were  taken  in  an  area  of 


Fig.  8.  Comparison  of  temperature  profiles  taken  at  the  locations  de¬ 
scribed  in  Figure  7. 


Fig.  9.  Vertical  variations  of  heat  content  in  the  leads  described  in  Fig¬ 
ure  7.  Total  heat  content  in  the  upper  2  m  of  the  water  column  was  2.0  MJ  m'1 
at  Site  A,  6.5  MJ  m~’  at  Site  B,  and  8.6  MJ  m at  Site  C. 
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scattered  floes  5-15  m  in  diameter  with  open  water 
areas  only  3-4  m  across.  Ice  concentration  was 
much  lower  (perhaps  only  50%),  but  the  area  was 
15-20  km  away  from  open  water.  There  was  typi¬ 
cally  a  strong  density  gradient  in  the  upper  70-100 
cm  due  to  meltwater  input,  but  salinities  were  usu¬ 
ally  above  15VM,  suggesting  significant  vertical 
mixing.  Salinities  below  1  m  were  generally  a  few 
parts  per  thousand  lower  than  interior  values,  pre¬ 
sumably  reflecting  the  entrainment  of  surface 
meltwater.  Heat  content  in  the  upper  2  m  was  2-3 
times  larger  (depending  on  the  distance  from  an 
ice  wall)  than  at  the  interior  site.  The  reason  for 
the  greater  heat  content  was  probably  the  greater 
amount  of  open  water,  although  there  could  also 
be  differences  in  the  incident  shortwave  radiation. 
Three  profiles  were  taken  on  22  July  across  a 
50-m-wide  lead.  Heat  content  was  substantially 
larger  than  in  the  small  leads  sampled  the  previous 
day  and  about  4  times  as  large  as  the  interior  site, 
indicating  that  both  lead  size  and  ice  concentra¬ 
tion  affect  the  heat  content.  There  was  also 
evidence  of  more  rapid  vertical  mixing,  probably 
due  either  to  windier  conditions  or  to  greater  rela¬ 
tive  motion  between  the  floes.  While  heat  content 
in  the  upper  50  cm  was  roughly  constant  across  the 
lead,  there  was  about  25%  less  heat  in  the  water 
column  on  the  downwind  side  of  the  lead.  Slightly 
lower  salinity  in  the  water  column  on  the  down¬ 
wind  side  suggests  that  the  reason  for  the  lower 
heat  content  was  more  efficient  transfer  of  heat  to 
the  floe  edge,  which  accelerated  the  melting. 


Radiation  measurements 

Incoming  shortwave  and  longwave  radiation 
were  sampled  at  Tent  Floe  every  10  minutes 
throughout  the  drift  (Figures  10  and  11).  Hourly 
totals  are  presented  in  Tables  5  and  6.  Shortwave 
fluxes  clearly  show  the  effects  of  variations  in 
cloudiness  superimposed  on  the  diurnal  fluctua¬ 
tions.  On  clear  days  such  as  1  and  3  July,  peak 
shortwave  fluxes  were  about  70%  larger  than  on 
cloudy  days.  Figure  10  indicates  that  cloud  cover 
tended  to  be  thinner  at  night,  even  when  there  was 
a  heavy  overcast  during  the  daytime.  The  diurnal 
signal  in  the  longwave  radiation  was  weak  and  the 
variability  seen  in  Figure  1 1  is  primarily  the  result 
of  changing  cloud  conditions.  The  minima  in  the 
curve  correspond  to  clear  sky  periods.  Because  of 
the  large  amount  of  water  vapor  in  the  air,  even  in 
the  absence  of  clouds,  fluctuations  in  incident 
longwave  radiation  were  only  on  the  order  of 
30%,  considerably  smaller  than  the  cloud  effects 
evident  on  incident  shortwave  radiation. 

Incident  spectral  irradiances  were  measured  pe¬ 
riodically  for  different  degrees  of  cloudiness.  Ex¬ 
amples  of  spectra  obtained  under  clear  and  cloudy 
skies  near  local  noon  are  shown  in  Figure  12.  Max¬ 
ima  and  minima  seen  in  the  two  curves  correspond 
to  minima  and  maxima  in  the  absorption  spectra 
of  atmospheric  water  vapor.  The  reduction  in  in¬ 
coming  solar  radiation  by  typical  MIZ  clouds  was 
roughly  a  factor  of  0.6  in  the  visible  and  up  to  a 
factor  of  3  in  the  infrared.  In  addition  to  observa¬ 
tions  at  the  drifting  station,  spectral  data  were 


Fig.  10.  Incoming  shortwave  radiation  over  Tent  Floe  during 
the  drift  phase. 
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Tabic  5.  Hourly  totals  of  Incoming  shortwave  radiation  (W/m1)  meas¬ 
ured  over  Tent  Floe  during  the  drift  phase. 

June  July 

Hour  (NST)  30  12345678 


Table  6.  Hourly  totals  of  incoming  longwave  radiation  (W/m1)  measured 
over  Tent  Floe  during  the  drift  phase. 


June  July 


Hour  (NST) 

30 

1 

2 

3 

4 

5 

6 

7 

8 

0100 

313 

313 

313 

315 

311 

313 

317 

309 

0200 

313 

315 

311 

315 

309 

313 

315 

309 

0300 

315 

317 

311 

315 

311 

313 

319 

311 

0400 

313 

319 

311 

315 

309 

313 

319 

313 

0500 

309 

322 

311 

306 

311 

317 

319 

313 

0600 

306 

319 

309 

311 

311 

317 

319 

309 

0700 

287 

322 

304 

304 

311 

317 

319 

309 

0800 

295 

319 

304 

309 

313 

317 

319 

0900 

287 

322 

269 

311 

311 

322 

324 

1000 

276 

324 

260 

309 

315 

319 

324 

1100 

276 

322 

263 

313 

313 

322 

322 

1200 

274 

319 

261 

313 

317 

324 

322 

1300 

274 

319 

261 

309 

322 

322 

322 

1400 

258 

319 

267 

311 

326 

324 

319 

1500 

254 

311 

267 

306 

322 

324 

317 

1600 

317 

249 

306 

289 

304 

322 

317 

315 

1700 

317 

247 

304 

306 

293 

322 

315 

313 

1800 

317 

245 

302 

317 

265 

326 

319 

293 

1900 

315 

258 

304 

317 

269 

324 

324 

289 

2000 

317 

282 

302 

315 

306 

319 

322 

293 

2100 

315 

282 

282 

317 

306 

315 

317 

306 

2200 

315 

293 

287 

317 

309 

313 

317 

306 

2300 

315 

306 

287 

315 

311 

311 

315 

306 

2400 

313 

311 

309 

315 

311 

313 

317 

309 

gathered  onboard  the  Polarstern  during  its  cruise 
to  the  Greenland  coast.  The  data  showed  an  even 
greater  reduction  in  incoming  solar  radiation 
under  cloudy  skies,  which  appeared  to  be  the  re¬ 
sult  of  lower  surface  albedos  and  reduced  back- 
scattering. 


Spectral  and  total  albedos  were  measured  over 
all  available  surface  types:  snow-covered  ice,  thick 
and  thin  first-year  ice,  multiyear  ice,  and  ponded 
ice.  Examples  for  snow  and  thin  (75  cm),  melting 
first-year  ice  are  shown  in  Figure  13.  The  upper 
curve  in  Figure  1 3  is  very  similar  to  previous  data 
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Fig.  12.  Incoming  spectral  irradiance  under  clear  sky  conditions 
(solid  curve)  on  3  July  and  under  cloudy  conditions  ( dashed 
curve)  on  7  July. 
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Fig.  13.  Spectral  albedos  observed  near  the  drift 
station.  The  solid  curve  is  for  a  10-cm-thick  snow  cov¬ 
er  over  first-year  ice  under  clear  skies  on  3  July;  the 
dashed  curve  is  for  melting  first-year  ice,  75  cm  in 
thickness,  under  cloudy  skies  on  7  July. 


we  acquired  in  the  Pt.  Barrow  area  over  snow-cov¬ 
ered  ice  in  the  initial  stages  of  melting.  The  depres¬ 
sion  of  the  albedo  below  0.9  in  the  500-nm  region 
suggests  the  presence  of  some  contamination  such 
as  soot  or  dirt  in  the  snow  layer.  The  values  for  the 
thin  first-year  ice  do  not,  however,  resemble  any 
of  the  cases  we  have  previously  studied.  The  de¬ 
crease  in  albedo  with  wavelength  is  very  similar  to 
that  of  cold,  bare  ice,  but  the  values  are  about  0. 1 
lower.  Albedos  measured  over  the  more  common 
blue  ice  are  comparable  with  the  thin  ice  values  in 
the  400-500-nm  region,  but  are  much  lower  in  the 
near  infrared.  It  is  believed  that  this  thin  ice  was 
very  young  (2-3  months)  and  that  the  difference  in 
optical  properties  was  related  to  higher  salinity  in 
the  ice. 
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Ocean  Currents  and  Temperatures  in  the 
Center  of  Fram  Strait  During  MIZEX  83 

KENNETH  HUNKINS 

Lamont-Doherty  Geological  Observatory  of  Columbia  University 
Palisades,  New  York  10964 


An  array  of  moored  oceanographic  recording 
instruments  was  placed  in  Fram  Strait  at  the  be¬ 
ginning  of  the  Marginal  Ice  Zone  Experiment  in 
1983  (MIZEX  83)  and  recovered  at  the  end  of  the 
experiment.  The  array  was  situated  near  the  center 
of  the  Strait  on  the  sill  separating  the  Arctic  Ocean 
from  the  Greenland  Sea  (Figure  1).  Three  moor¬ 
ings  were  arranged  in  a  triangle  with  the  fourth 
mooring  at  the  center  of  the  triangle.  Nominal 
spacing  was  10  km,  a  distance  chosen  as  close  to 
the  Rossby  radius  of  deformation  in  the  region. 
The  array  was  in  open  water  east  of  the  ice  edge. 

The  objective  was  to  monitor  mesoscale  oceano¬ 
graphic  fluctuations  in  open  water  near  the  ice 
edge.  An  oceanographic  front  is  generally  associ¬ 
ated  with  the  ice  edge,  and  two  extreme  hypothe¬ 
ses  might  be  postulated  about  mesoscale  motion  in 
the  area.  Such  motions  could  be  due  to  a  spectral 
cascade  of  energy  generated  in  boundary  layers  by 
the  mean  flow.  In  this  case,  spectra  in  the  time  do¬ 
main  would  be  expected  to  show  a  monotonic  de¬ 
crease  in  energy  from  low  to  high  frequencies,  rep¬ 
resenting  a  cascade  of  geostrophic  turbulence.  An 
alternate  hypothesis  would  be  that  instability 


processes  in  the  mean  shear  across  the  front  favor 
the  growth  of  certain  discrete  frequencies  and 
wavenumbers.  In  this  case  the  recordings  would  be 
expected  to  show  one  or  more  spectral  peaks  in  the 
subinertial  frequency  range. 

Four  subsurface  moorings  were  deployed  in 
water  depths  of  about  2400  m  on  16-17  June, 
1983,  from  M/V  Polarbjdrn  in  an  array  nominally 
centered  at  78  °30  'N,  2  °E  (Table  1).  The  array  lo¬ 
cation  was  south  of  the  Molloy  Deep,  where 
depths  exceed  4000  m,  and  north  of  a  northwest- 
southeast  trending  ridge  that  rises  to  a  depth  of 
about  1700  m.  Bottom  slopes  were  gentle  in  the 
area  of  the  array.  All  of  the  moorings  and  instru¬ 
ments  were  retrieved  41  days  later  on  27  July, 
1983.  One  of  the  moorings  was  then  redeployed  to 
record  over  the  next  year.  It  will  be  retrieved  dur¬ 
ing  MIZEX  84.  Aanderaa  RCM-5  current  meters 
were  used  to  record  current  speed  and  direction,  as 
well  as  temperature  and  depth.  Manufacturer’s 
specifications  for  these  CTD  instruments  are: 
Temperature  resolution  ±0.01  °C 

Current  speed  accuracy  ±2%  of  actual  speed 

Current  direction  resolution  ±0.35° 


Table  1.  Mooring  locations  and  depths. 


Fig.  1.  Location  of  MIZEX  83  mooring  array.  (Bathymetry  from  Johannessen  and  Johannessen,  1983.) 
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All  instruments  were  calibrated  before  the  MIZEX 
83  cruise.  Although  the  instruments  were  equipped 
with  conductivity  sensors,  range  settings  were  in¬ 
correct  for  this  region  and  no  salinity  data  were 
obtained.  Current  meters  were  suspended  at  nomi¬ 
nal  depths  of  30  and  130  m  below  the  surface  on 
all  four  moorings.  Additional  instruments  were 
suspended  at  depths  of  1000  m  below  the  surface 
and  100  m  above  the  bottom.  Data  were  sampled 
every  10  minutes. 

The  magnetic  recording  tape  failed  to  transport 
on  four  of  the  twelve  instruments  deployed.  The 
malfunction  was  later  found  to  be  caused  by  ex¬ 
cessive  use  of  tape  head  cleaner,  which  did  not  en¬ 
tirely  evaporate  before  the  tape  was  mounted. 
This  resulted  in  adhesion  of  tape  to  the  tape  head. 
Another  instrument  was  accidentally  fouled  in  the 
mooring  cable  during  launching  so  that  it  was  un¬ 


able  to  swivel.  Temperature  data  are  available  for 
that  instrument,  but  current  data  are  invalid.  On 
three  other  instruments  digital  noise,  apparently 
due  to  poor  functioning  of  the  digitization  mech¬ 
anism  in  these  three  meters,  reduced  the  length  of 
useable  record. 

The  moorings  were  situated  in  open  water  30 
km  east  of  the  ice  edge,  which  trends  nearly  north 
and  south  in  this  region.  An  oceanographic  profile 
was  made  from  M/V  Polarbjdrn  by  the  University 
of  Bergen  Oceanographic  Group,  with  a  Neil- 
Brown  CTD  instrument,  just  before  the  moorings 
were  retrieved  (Figure  2).  Sigma-r  was  about  27.2 
at  the  surface  and  increased  rapidly  with  depth  in 
the  upper  100  m  so  that  at  130  m  it  was  about  27.8. 
Below  100  m  it  increased  slowly  and  was  still  less 
than  28.0  at  2000  m.  There  was  a  thin  layer  of 
warm  surface  water  at  about  5°C  that  extended  to 
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Table  2.  Mean  values  and  standard  deviations  for  tem¬ 
perature  and  current. 


Mooring 

number 

Instrument 

depth 

(m) 

Record 

length 

(days) 

T 

(°C) 

U 

(cm/s) 

V 

(cm/s) 

§ 

1 

52 

26 

3.62 

±0.56 

-1.7 

±8.2 

0.4 

±9.4 

• 

151 

29 

1.6 

±5.3 

-5.0 

±7.3 

2 

61 

41 

2.51 

±0.36 

- 

- 

• 

1015 

41 

-0.77 

±0.08 

1.1 

±2.5 

-4.4 

±4.2 

4 

55 

41 

2.56 

±0.63 

-2.0 

±6.6 

-7.5 

±8.2 

158 

30 

1.62 

±0.40 

-2.0 

±5.8 

-6.5 

±7.3 

• 

1014 

41 

-0.74 

±0.08 

0.5 

±2.4 

-5.2 

±3.8 

- 

2296 

41 

-1.00 

±0.04 

1.2 

±3.0 

-7.5 

±3.9 

__ 

rectly  with  the  other  two  means  of  2.51  and 
2.56  °C,  which  are  satisfyingly  close  to  each  other 
in  value.  Their  large  standard  deviations,  ±0.36 
and  ±0.63  °C,  reflect  the  fact  that  these  instru¬ 
ments  are  located  in  a  thermocline,  so  that  any 
vertical  water  motions  produce  large  recorded 
temperature  variations. 

The  two  instruments  at  the  150-m  level  recorded 
mean  temperatures  of  1 .67  and  1 .62  °C,  again  rela¬ 
tively  close  in  value,  as  are  also  the  two  means  of 
-0.77  and  -0.74 °C  at  the  1000-m  level.  There  was 
only  one  record  from  100  m  above  the  bottom  and 
the  mean  at  that  level  was  -1.0°C,  quite  close  to 
the  temperature  at  2000  m  observed  for  the  CTD 
profile.  The  standard  deviation  of  ±0.04°C  near 
the  bottom  indicates  the  constancy  of  temperature 
at  that  level  and  also  gives  some  measure  of  the  re¬ 
peatability  of  the  temperature  measurements. 

Mean  curents  at  four  levels  in  the  water  column 
from  near  the  surface  to  near  the  bottom  all  have  a 


southward  component  (Figure  3).  Mean  speeds  are 
relatively  high,  ranging  from  5.2  to  7.6  cm/s.  Al¬ 
though  there  are  shears  between  the  different 
levels,  it  is  perhaps  the  agreement  rather  than  the 
differences  between  these  means  that  should  be 
emphasized.  The  general  similarity  in  speed  and 
direction  between  the  different  levels  suggests  that 
the  mean  flow  has  a  large  barotropic  component. 
The  time  series  do  show  some  long-period  trends 
extending  over  several  weeks,  and  means  over 
longer  periods  might  have  shown  different  results. 
The  results  do  indicate  that  there  must  have  been  a 
large  northward  flow  in  another  part  of  the  Strait 
during  this  period  to  compensate  for  the  large  ob¬ 
served  southward  transport. 

The  MIZEX  83  mooring  results  represent  the 
combined  efforts  of  many  people.  The  moorings 
were  deployed  under  the  supervision  of  Ivars 
Bitte.  Jay  Ardai,  Barry  Allen,  and  Peter  Bruch- 
hausen  assisted  with  the  launching.  Recovery  was 


Fig.  3.  Mean  current  vectors  at  Mooring  4.  Depths  in  meters.  41 -day  means,  except  for  the  158-m  record,  which 
30-day  mean. 
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Our  motive  for  recording  whitecap  coverage  of 
the  sea  at  high  latitudes,  and  particularly  in  the  vi¬ 
cinity  of  the  ice  pack  during  MIZEX  83,  was  to  ex¬ 
tend  the  sea- water  temperature  ( Tw )  range  repre¬ 
sented  in  our  whitecap  (H7),  windspeed  (U)  data 
file  so  that  we  could,  by  application  of  a  suitable 
statistical  methodology  (e.g.  O’Muircheartaigh 
and  Monahan,  1983a),  obtain  an  improved 
W(U,  T.)  expression.  Our  hypothesis  when  we  un¬ 
dertook  this  study  to  determine  the  influence  on 
W  of  Tw,  as  distinct  from  the  influence  of  atmo¬ 
spheric  thermal  stability  or  the  stability-related 
measure  T.  -  T.  (where  Ta  is  the  air  temperature), 
was  that  W  would  decrease  with  increasing  T.. 
This  hypothesis  was  based  on  the  following  as¬ 
sumptions: 

1)  For  a  given  wind  speed  (and  stability), 
that  a  uniform  number  of  waves  would 
break  per  unit  area  of  sea  surface  per  unit 
time,  regardless  of  Tw; 

2)  That,  for  any  Tw,  the  average  breaking 
wave  would  entrain  the  same  amount  of  air 
to  the  same  depth;  and 

3) That  the  bubble  spectrum  resulting  from 
wave  breaking  would  be  independent  of  7V 
When  assumptions  2  and  3  are  combined  with 

Blanchard’s  early  observations  (personal  com¬ 
munication),  which  showed  that  for  any  given 
bubble  size  the  rise  velocity  of  the  bubble  increases 
significantly  when  there  is  a  modest  increase  in  Tw, 
the  conclusion  is  reached  that  cold-water  white- 
caps  persist  longer  than  warm-water  ones,  since 
the  whitccaps  in  cold  water  are  replenished  by 
bubbles  from  below  for  a  longer  period  than  are 
whitecaps  in  warm  water.  This  conclusion,  when 
taken  together  with  assumptions  1  and  2,  leads  di¬ 
rectly  to  the  aforementioned  hypothesis. 

We  now  have  at  hand  the  first,  preliminary  re¬ 
sults  from  the  analysis  of  some  of  the  MIZEX  83 
whitecap  photographs.  While  the  preliminary  re¬ 
sults  shown  in  Figure  1  are  based  on  the  work  of 
only  one  analyst  and  ignore  the  influence  of  stabil¬ 


ity  and  fetch,  they  have  already,  when  considered 
in  conjunction  with  some  of  our  recent  whitecap 
simulation  tank  findings,  caused  us  to  revise  our 
working  hypothesis  drastically.  When  the  points 
on  Figure  1  are  compared  with  the  solid  W(U) 
curve  (Monahan  and  O’Muircheartaigh,  1980) 
based  on  warm  ocean  observations  reported  in 
Toba  and  Chaen  (1973)  and  Monahan  (1971),  or 
with  the  dashed  curve  (O’Muircheartaigh  and 
Monahan,  1983b,  Eq.  9)  based  on  the  relatively 
cool-water  (12.5°-I4°C)  JASIN  whitecap  observa¬ 
tions  (Monahan  st  al.,  1981),  the  suggestion  arises 
that  for  a  given  U,  W  is  less  for  cold  seas  than  for 
warmer  waters.  Indeed,  these  preliminary  MIZEX 
83  low  W  results  are  reminiscent  of  the  fresh-water 
whitecap  observations  reported  by  Monahan 
(1969),  even  to  the  relatively  high  wind  speed  for 
the  onset  of  whitecapping — the  Beaufort  velocity 
discussed  in  O’Muircheartaigh  and  Monahan 
(1983b). 

While  these  surprising  initial  MIZEX  results 
were  coming  to  light,  some  equally  intriguing  re¬ 
sults  were  being  obtained  with  the  U.C.G.  whitecap 
simulation  tank.  The  amount  of  space  charge  pro¬ 
duced  during  the  decay  of  the  whitecap  formed  by 
a  ‘standard’  breaking  wave  in  the  tank  was  found 
to  increase  several  fold  as  the  temperature  of  the 
sea  water  in  the  tank  was  increased  from  5°  to 
25  °C  (see  Figure  5,  Bowyer,  1983).  This  finding 
lent  itself  to  one  of  two  explanations:  either  the 
space  charge  per  aerosol  droplet  increases  marked¬ 
ly  with  increases  in  Tw,  or  the  number  of  droplets 
(particularly  jet  droplets  in  this  context)  produced 
per  whitecap  increases  significantly  with  elevation 
of  T..  Subsequent  experiments,  in  which  a  Royco 
Model  225/241  aerosol  particle  counter  was  used 
to  measure  the  number  of  droplets  injected  into 
the  hood  of  the  tank  as  a  result  of  ‘standard’ 
breaking  waves  produced  at  various  7V s,  showed 
that  the  actual  number  of  jet  droplets  produced 
per  breaking  wave  increases  markedly  with  in¬ 
creasing  Tw.  While  it  has  not  yet  been  confirmed 


Fig.  1.  The  fraction  of  the  sea  surface  covered  by  whitecaps,  W,  vs 
the  wind  speed  at  10  m  elevation,  U;  preliminary  results  of  MIZEX  83 
observations.  Open  squares  12.5°C  <  T.  s  14°C;  filled  triangles 
S.OeC  <  T.  s  12. 5°C;  filled  circles  T»  £  5.0°C.  Each  point  is  based 
on  the  analysis  of  8  to  22  photographs.  This  preliminary  figure  based 
on  the  analysis  of 346  photographs  out  of  some  1040  taken  from  F/S 
Polarstern.  The  solid  and  dashed  curves  are  described  in  the  text. 
N.B.  Points  near  bottom  of  panel  with  downward  pointing  arrows 
have  W  values  less  than  0.0001.  Many  correspond  to  intervals  with 
no  whitecaps. 


experimentally,  it  is  probable  that  this  increase  in 
jet  droplet  production  with  increasing  Tw  is  due  to 
the  fact  that,  as  Tw  increases,  the  number  of  bub¬ 
bles  increases  and  the  bubble  size  spectrum  peak 
shifts  to  smaller  bubble  radius,  probably  as  a  re¬ 
sult  of  reduced  viscosity  (smaller  bubbles  are 
much  more  efficient  generators  of  jet  droplets,  as 
was  shown  by  Blanchard,  1963). 

If  this  interpretation  of  the  tank  results  is  cor¬ 
rect,  then  assumption  3  (and  perhaps  also  assump¬ 
tion  2)  listed  above  is  wrong.  We  thus  are  left  with 
a  new  working  hypothesis:  that  cold-water  white- 
caps  persist  for  briefer  periods  of  time  than  warm- 
water  whitecaps  because,  in  spite  of  the  greater 
viscosity  of  the  cold  water,  the  larger  bubbles  that 
replenish  the  cold-water  whitecaps  reach  the  sur¬ 
face  sooner  than  the  small  bubbles  that  sustain 
warm-water  whitecaps.  A  concomitant  of  this  hy¬ 
pothesis  is  that  for  the  same  U,  W  is  less  in  cold 
seas  than  in  warm.  Indeed,  this  argument  is  essen¬ 
tially  the  same,  with  regard  to  the  effect  of  bubble 
spectra  on  W(U ),  as  that  put  forth  to  explain  the 
relatively  low  W  values  observed  on  the  North 
American  Great  Lakes  as  compared  to  those  ob¬ 
tained  from  oceanic  observations  (Monahan  and 
Zietlow,  1969;  Monahan,  1971). 

We  hope  shortly  to  use  a  video  area  analyzer  in 
conjunction  with  our  MIZEX  83  whitecap  tapes  to 
test  this  new  hypothesis  as  regards  the  short  life¬ 
times  of  cold-sea  whitecaps.  We  look  forward  to 
collecting  additional  video  and  photographic  rec¬ 
ords  during  MIZEX  84  that  will  enable  us  to  refine 
the  W(U,Tw)  expression  that  we  will  initially  base 
on  the  MIZEX  83,  STREX  (Doyle,  1984),  JASIN, 
and  BOMEX  whitecap  observations.  The  resulting 
W(lI,Tw)  expression  can  then  be  introduced  into 
our  sea  surface  aerosol  generation  model  (Mona¬ 
han  et  al.,  1982,  1983a),  and  the  improved  model 
can  be  used  to  obtain  revised  estimates  of  global 
sea-to-air  salt  fluxes. 

The  research  discussed  in  this  note  is  sponsored 
by  the  U.S.  Office  of  Naval  Research  via  grant 
N00014-78-G-0052.  The  MIZEX  83  whitecap  ob¬ 
servations  were  recorded  aboard  the  F/S  Polar- 
stern  of  the  Alfred  Wegener  Institut,  Bremerha- 
ven.  The  assistance  of  Prof.  E.  Augstein,  the  other 
scientists,  and  the  crew  of  the  Polarstern  is  grate¬ 
fully  acknowledged.  The  whitecap  results  present¬ 
ed  in  Figure  1  are  based  on  the  photographic 
analyses  of  Miss  M.  Higgins.  The  extrapolation  of 
the  winds  to  10-m  elevation  values  was  done  by 
Dr.  M.C.  Spillane. 
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